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ABSTRACT 


This  report  details  the  final  test  results  of  coioputer  program 
WARP.  This  program  generates  coefficients  of  a  least-squares 
LORAN  grid  warpage  model  for  use  in  the  AN/ARN-101  Digital 
Avionics  System.  The  coefficients  enhance  the  capability  of  the 
AN/ARN-101  to  accurately  convert  from  LORAN  time  differences  to 
Latitude-Longitude  pairs.  Test  results  include,  input  error 
cliecks,  boundary  jump  analysis,  altitude  effects  and  in))ut  data 
de  ns i ty . 


FIUJUj 


TABLE  OF  CONTENTS 


PAGE 


1.0  INTRODUCTION  1 

1.1  Purpose  1 

1.2  Scope  1 

1.3  Test  Overv i ew  1 

1.4  General  Test  Procedures  3 

2.0  APPLICABLE  DOCUMENTS  5 


2.1  MethodofTest  5 

2.2  User ' s  Manual  5 

2.3  AN/ARN-101  Specifications  5 

3.0  TEST  OBJECTIVES  6 


3.1  OBJECTIVE  1  (Bad  Input  Data) 

3.1.1  Description 

3.1.2  Test  Procedure 

3. 1.2.1  Program  Modification 

3. 1.2.2  Data  Requirements 

2.1.3  Test  Results 

3.1.4  Evaluation  Criteria 

3.1.5  Test  Evaluation 

3.1.6  Recommendations 

3.2  OBJECTIVE  2  (Program  Constants) 

3.2.1  Description 

3.2.2  Test  Procedure 

3.2.2. 1  Program  Modification 

3. 2. 2. 2  Data  Requirements 

3.2.3  Test  Results 

3.2.4  Evaluation  Criteria 

3.2.5  Test  Evaluation 

3.2.6  Recommendations 

3.3  OBJECTIVE  3 

3.3.1  Description 

3.3.2  Test  Procedure 


6 

6 

6 

6 

6 

6 

6 

6 

7 
3 
3 

8 
8 
8 
8 

12 

12 

12 

13 

13 

13 


V 


PAGE 


3.3.2. 

L 

Program  Mooi f i cat i on 

14 

3.3.2. 

?. 

Data  Reqni rements 

14 

3.3.3 

Test  Resul ts 

14 

3 .3.4 

Tva 1 uat  i  on  Criteria 

14 

3 . 3  .  b 

Test  Eva  1 uat i on 

14 

3 . 3 .  r, 

Recoitunendat  ions 

17 

3.4 

OBJECTIVE  4 

(Input  Data  Oensity/Size  of  PCA) 

18 

3.4.1 

Description 

18 

3.4.2 

Test  Procedure 

18 

3.4.2. 

1 

Program  Modification 

18 

3.4.2. 

Data  Pequ i rements 

18 

3.4.3 

Test  Results 

18 

3.4.4 

Evaluation  Criteria 

18 

3.4.3 

Test  Evaluation 

18 

3.4.  n 

Rec  oiiime  nda  t  i  o  ns 

19 

3.') 

OBJECTIVE  5 

(TD  Predictions  in  Void  Areas) 

20 

3.3.1 

Description 

20 

3.3.2 

Test  Procedure 

20 

3.3.:^. 

[ 

Program  Modification 

20 

3.32. 

V 

Data  Requi  reuients 

20 

3.5.3 

lest  Results 

20 

3.3.4 

Evaluation  Criteria 

20 

3.3.3 

Tost  Evaluation 

20 

3.3.3 

Rocommo  nda  t i ons 

20 

3.3 

OBJECTIVE  6 

(Altitude  TD  Predictions) 

21 

3.6.1 

Desert  jjtion 

21 

3  ,  3  . 

Test  Procedure 

21 

3.3.2. 

1 

Program  M  o  d i f i c  a  t i 0  n 

21 

3.3.2. 

y 

Da  1,  a  R eq  u  i  r erne n t  s 

21 

'  .  3  .  1 

Tost  Results 

21 

i .  3 . 4 

Evaluation  Criteria 

23 

3.3.3 

Test  Evaluation 

23 

3.3.  3 

Rec  omine  nda  t  i  on s 

24 

3.  1 

OBJECTIVE  7  (Boundary  Conditions) 

33 

3  .  / .  1 

Dose  r i j  t i 0  n 

33 

3  .  /  .  2 

Test  P  r  0  c  0  n  u  r  e 

33 

3.7,2. 

1 

L 

Prpyrai.t  Modification 

33 

3.7.'?. 

2 

Data  Req  u  i  roine n t s 

33 

1.1.1 

Test  Results 

33 

3.7.4 

Evaluation  Criteria 

37 

3.7,3 

Test  Evaliiation 

37 

1.1  .h 

Recoiiii.ienda  t  i  ons 

37 

3  .  . 

OBJECTIVE  H  (Average  and 

Saltwater  Impedances ) 

39 

3.3.  1 

Dose  r i pt i on 

39 

3.3.  ' 

Test  Procedure 

39 

3  ,  .'1 .  2  . 

1 

P  r  0  y  r  a  1.!  Modification 

40 

VI 


PAGE 


3. 8. 2. 2  Data  Requirements  40 

3.8.3  TestResults  40 

3.3.4  EvaluationCriteria  40 

3.8.5  Test  Evaluation  40 

3.8.6  Recommendations  40 

3.9  OBJECTIVE  9  (Photo  Reconnaissance  44 

Data  Validation) 

3.9.1  Descri pt i on  44 

3.9.2  Test  Procedure  44 

3.9.2. 1  Program  Modification  44 

3. 9. 2. 2  OataReqiurenebts  44 

3.9.3  TestResults  44 

3.9.4  Evaluation  Criteria  44 

3.9.5  Test  Evaluation  44 

3.9.6  Recommendations  44 

3.10  OBJECTIVE  10  (Software  Program  45 

Document  at  ion) 

3.10.1  Descri pt i on  45 

3.10.2  TestProcedure  45 

3.10.2.1  Program  Modification  45 

3.10.2.2  Data  Requirements  45 

3.10.3  TestResults  45 

3.10.4  Evaluation  Criteria  45 

3.10.5  TestEvaluation  45 

3.10.6  Recommenda t i ons  45 

4.0  CONCLUSIONS  46 

5.0  RECOMMENDATIONS  48 

Appendix  A  to  ran  Warp  age  Data  Base  Requirements  A-1 

Appendix  B  Statement  of  Work  B-1 


LIST  OF  !A3LES 


lABLES 


TITLE 


PAGE 


I  WARP  Input  Parameters  4 

II  Constant  Ranges  and  Nominal  Values  8 

III  Mean  TO  Errors  as  Functions  of  Program  9 

Constants 

IV  Mean  TO  Errors  as  Function  of  Vertical  Lapse  9 

V  Reference  Spheroids  13 


VI 


VI  I 


VI  It 


IX 


X 
XI 
XI  [ 

x:  1 1 

XIV 

XV 


Positional  Error  Oue  to  Inconsistent  15 

Spheroid  Selection  in  the  Warpage  Coefficient 
Generation  Process 

Positional  Error  Due  to  Incorrect  15 

Spheroid  Selection  in  the  AM/ARN-101 

Loran  Warpage  Model  Accuracy  as  a  19 

Function  of  Point  Density 

LORAN  Warpage  Positional  Error  21-22 

Area  2  Boundary  Conditions  35 

Area  3  Boundary  Conditions  35 

Area  4  Boundary  Conditions  36 

Area  5  Boundary  Conditions  36 

Constant  Impedance  Values  39 

Position  Error  (in  feet)  of  Impedance  Models  42 


LIST  OF  FIGURES 

FIGURES  TITLE  PAGE 


1  F-4E  Aircraft  Oata  Collection  Flight  Profile  2 

2  Mean  TO  Errors  as  a  Function  of  Vertical  10 

Lapse 

3  Mean  Error  {Feet)  as  a  Function  of  Vortical  11 

Lapse 

4  Spheroid  Probaoility  Functions  16 

5  Warpage  Contour  Plot  (usee)  -  Slave  A,  IK  PS 

6  Warpage  Contour  Plot  (usee)  -  Slave  A,  5K  26 

7  Warpage  Contour  Plot  (usee)  -  Slave  A,  lOK  27 

8  Warpage  Contour  Plot  (usee)  -  Slave  A,  15K  28 

9  Warpage  Conto'.r  Plot  (usec)  -  Slave  8,  IK  20 

10  Warpage  Contour  ^lot  (usec)  -  Slave  B,  5K  30 

11  Warpage  Contour  Plot  (usec)  -  Slave  3,  lOK  31 

12  Warpage  Contour  Plot  (usec)  -  Slave  3,  15K  32 

13  Boundary  Junip  Errors  34 

14  Constant  Impedance  Probability  '^unctions  at  IK  43 


1.9  INTRODUCTION.  This  document  is  the  Test  Report  for  the  Warpage  Coefficient 
Ceneration  Program  (WARP).  This  off-line  computer  support  program  uses  ri3al  world 
Loran  time  difference  (TD)  measurenents  to  produce  i.oran  grid  warpage  correction 
factors  to  enable  precise  AN/ARN-101  navigation  and  v;eapon  delivery. 

1.1  PURPOSE.  The  test  of  the  Warpage  Coefficient  Generation  Program  was  designed 
to  define  data  requirements,  evaluate  program  utility  and  measure  toran 
navigational  accuracy.  The  test  results  provide  a  basis  for  fomiulutinj 
guidelines  for  the  operational  use  of  the  coefficient  generation  progran  in  a 
tactical  environment. 

1.?  SCOPE.  The  extensive  data  base  collected  during  Afl/ARN-lOl  flig!\t  testing  at 
Cglin  AF3,  Florida  was  used  to  exercise  WARP.  The  prograiii  input  dati  consists  of 
four  altitude  dependent  sets  of  precise  geodetic  latitude  longitude  (!,Vr/!_()N) 
coordinates  and  the  corresponding  raw  Loran  time  difference  measure nents  within 
the  Southeast  USA  Loran-C  Chain  coverage  a'-ea.  "he  test  defines  Loran  accuracy  as 
a  function  of  WARP  parameters,  input  data  density,  altitude,  and  size  of  the  j)rime 
coverage  area  (PCA).  The  primary  goal  of  this  test  activity  is  to  quali^'y  input 
data  requirements  and  evaluate  Loran  navigational  accuracy.  The  following  test 
objectives  are  addressed: 

a.  Identification  of  bad  input  data 

b.  Selection  of  program  constants 

c.  Improper  spheroid  selection 

d.  Effects  of  input  data  density/size  of  PCA 

e.  TO  predictions  in  void  areas 

f.  Altitude  TO  predictions 

g.  Boundary  cofiditions 

h.  Average  and  saltwater  impedances 

i.  Photo  reconnaissance  data  validation 

j.  Software  program  documentation 

1.3  TEST  OVERVIEW.  The  data  base  fo^  this  test  effort  was  established  in 
‘fulfilment  of  Objective  IB  of  the  Annament  Division  Test  Directive  No.  PB'VirAO,?, 
AN/ARN-101  i^ave  Tack  Interoperability  Flight  Test  (F-4  Advanced  Avionics  P'naso 
II).  This  objective  addressed  the  mapping  of  the  Southeast  USA  Loran-C  Chain  for 
the  purpose  of  generating  Loran  warpage  coefficients.  Four  AN/ARN-101  equipped 
F-4E  airborne  missions  were  used  to  collect  data  at  altitudes  of  I'iOO  ft  MSL ,  5000 
ft  MSL,  10000  ft  MSL,  and  15000  ft  MSL.  Mosaic  flight  profiles  were  flown 
throughout  the  PCA  as  shown  'n  Figure  1. 


'"'■'-.L ,  (I'jous  du'jl  "'’S-ib  radar  and  se''ec';ive  ;)'nol,othaodol  i  I;.?  covorjj,'  I 

ior  accurate  s.)ace  positioning  of  the  airc-'aft.  fhe  aircraft  post-'^l  i  j  it 
’  nstruiiientdt  ion  data  with  nerged  dual  radar/ph  jtotheodol  i  ta  pos  ’  t  i  )ni  n-j  data  ./ea- 
t  ime  sequericed  and  merged  to  provide  the  required  WAP.'’  data,  hata  retuci  io'i 
techniques  provided  coordinate  paired  LAT-^ON's,  fOs,  and  altitude  da'  i  ,))int>  for 
each  of  the  four  flight  levels,  fhe  selection  of  data  jioints  fo''  a,  ■>,  'dti’" 
input  was  bounded  by  a  3-space  aircraft  position  error  of  loss  than  inn  t',.  inly 
selected  points,  as  necessary  to  satisfy  the  5  nn  di  s  tr  i  Iju  t  ’  on  des  '  r  3')  i  1  i  ty ,  were 
chosen  with  piosition  errors  of  103-200  feet. 

"^he  revised  'WARP  and  the  four  altitude  depc’iidont  d.ata  leases  are  c  irr.ettly 
•"osident  on  the  fglin  A,-3  tOC  hh'lC  cO'.iputer.  'I'o  each  data  file  is  attac'n.'d 
iddi'.'onal  WARP  parameters  which  ih.'fine  th-}  test  area  'ooundar  i  .?s ,  torin  cha  i 'i 
definitions,  the  earth  spiieroid  aodel  definition,  and  noaiin.il  'ragi-d  1  coristan'  .. 

A  listing  of  these  program  parameters  is  given  mi  ^al)le  I.  A  de!,a'"'i.'d  inaly'.i'. 
Loran  gr-'d  warpage  and  grid  warpage  correction,  are  areJ  by  ‘Ir  J.  ...  ''oward,  t.he 
Corporation,  si’ec  i '' i  cal  ly  acdresses  t'ne  Cglin  activity  and  is  'nclud'-'!  ,!■, 
Appendix  A. 

1.4  GENERAL  TEST  PROCEDURES.  A  statistical  evalu.ator  was  used  to  exercise  the 
v^arpage  inodel  against  a  'oascline  set  of  measur.ed  LA^-LON'/TD  coordinat'.'  pam-s.  "he 
evaluator  consisted  of  twa  baseline  data  sets  (A  and  ?<)  ^o''  eacli  o'"  the  requir.-,: 
test  altitudes.  LAT-LON/tD  and  altitude  c'ata  points  we'"e  selected  to  develop 
equ’ valent  sets  with  a  density  distribution  of  one  point  per  G  ma  ce'l  pr-jvi.dinj 
38  cell  points  within  the  Eglin  ?CA.  The  two  sets  -.vere  doternined  ^,0  be 
statistically  equiv.alent  with  respect  to  WAR'’  mean  tamar,  standarii  devi.a'Jon,  .and 
r.ange  conpusat  i  ons.  Once  i  nterchangeabl  e,  one  sia'  was  id'-ant  i  f  i  la  .1  as  '.li>  'WAR'' 

■n;!.:t  data,  and  trae  other  laecame  ‘.h.e  evaluation  set  used  to  nia  lOiin’inj 

a  c  c  u  a  cy . 
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2.0  APPLICABLE  DOCUMENTS. 


The  following  documents  are  api'Iicahle  as  defined. 

2.1  METHOD  OF  TEST. 

The  basis  for  this  report  is  the  Method  of  Test  (MOT)  for  the  Loran  Warpaye 
Coefficient  Generation  Program  (WARP).  The  MOT  identifies  ten  (10)  objectives  for 
testing  the  software  program.  It  was  prepared  by  the  Electronic  Systeins  Division 
Operating  Location  -AF  (ESD/OL-AF)  at  Eglin  AFB,  Florida.  The  MOT  is  dated  20  Doc 
79  and  is  available  through  ESD/OCN-1,  Hanscom  AFB,  MA. 

2.2  USER'S  MANUAL. 

The  Warpage  Coefficient  Generation  Program  was  developed  by  Lear  Sieyler, 

Inc.  (LSI)  and  provided  to  ESO.  The  program  source  listing  v;as  accompanied  by  a 
preliminary  User's  Manual  LSI  YV1007,  13  May  76.  The  User's  Manual  wa;  updated  by 
ESD/OCN-1.  At  the  time  of  testing  this  manual  (OCN  79-305,  1  Nov  79)  coiiiprisod 
the  most  current  program  documentation. 

Results  of  ESD  testing  include  a  revised  Fortan  program,  a  new  User's  Manual , 
and  an  Analyst's  Manual  (See  Objective  10). 

2.3  AN/ARN-101  SPECIFICATION. 

The  warpage  correction  reguirem»nts  are  determined  by  the  Computer  Prograiii 
Development  Specification  for  Operational  Flight  Program,  RF-AC  (F-4L)  Digital 
Modular  Avionics  Systeiii  AN/ARN-101  (V),  Contract  Number  F19623-76-C-002A ,  Oocuiiient 
Number  CB1001-004A  (CBlOOl-OlOA) . 


3.0  riST  OBJECTIVES. 


3.1.1  DESCRIPTION.  DemonsLrate  that  the  program  ic  rapable  of  identifying  bad 
iripiit  data. 

3.1.’  TEST  PROCEDURE.  Examine  WARP  to  dotenni ne  the  capability  to  detect  had  or 
erroneous  data. 

a.l.’.l  PROGRAM  MODIFICATION.  None 

3.’.’.’  DATA  REQUIREMENTS.  Nominal  program  input  data  sets. 

3.1.3  IlST  results.  Nominal  prograi.i  outputs. 


3.1.4  I. VALUATION  CRITERIA.  The  input  da: a  set  consists  of  correlated  Loran  time 
differences  and  geographical  latitude/lonjitudes.  The  definition  of  bad  or 
ert-oneous  data  is  limited  to  mean  those  data  points  where  the  ..leasured  time 
lii;erences  tit)  not  correspond  witi  their  associated  LAT/tON  for  reasons  ot'ner  tiian 
warpage.  Sinct!  Loran  war,. age  is  defined  lo  be  the  divergence  of  the  Loran 
refi'rt.'nce  relative  to  a  given  earth  reference  due  to  jihysical  eartli  properties,  it 
is  (jnife  difficult  ti)  soiiarate  wanidge  from  measurement  error. 

Within  WARP  the  only  error  detection  lechanisi.i  is  a  test  to  determine  if  a 
c'xii'uted  impedance  lies  within  the  real  world  range  of  .001056  Z  .08. 
Imped.ances  t'nat  fall  outside  this  range  are  not  found  in  nature  and  are  therefore 
considered  as  bad  data  points. 

3.1.5  ILST  EVALUATION.  The  computer  program  estimates  an  impedance  using  a 
■Jew' o'l-Raphson  iterative  sciie.ue.  When  an  interim  value  violates  the  proscribed 
limiLu,  tile  iteration  ceases  and  results  in  a  non-zero  error  residual.  This 
non-c mvorgence  error  residual  is  used  as  the  flag  to  identify  bad  data  points. 
Since  tfiere  are  three  unknown  impedances  and  only  two  time  difference  equations, 

!  seli’cted  a;)i)roacli  within  WARP  sots  the  master  station  impedance  to  a  constant 
value  and  iterates  the  slave  impedances  to  reach  a  simultaneous  equation  solution. 
A  po  r  choice  of  master  impedance  cun  cause  the  iteration  scheme  to  attempt  to 
'.ulr':'  ;ion-rea  1  izu.)I  e  i'.ijiedance  values  fc-  either  slave.  A  look-up  table  limit 
val'i'  m,  then  reac'n.ed  and  a  res.cual  TD  error  resalts,  thus  identifying  "bad” 
da'.a.  A  rjiige  of  master  impedance  valuer  can  generate  a  multiple  nuinber  of  valid 
s.)Iu'.  1  ons  tiidt  r.iay  or  ..lay  not  re|.resent  real  world  values.  Conversely,  a  ;)oor 
choir.'  M,ay  also  generate  TD  resicuals.  (See  also  Objective  2.). 

It  is  therefore  concluoed  thc.t  an  errar  residual  flag  is  not  necessarily  an 
indiiMtor  of  "bad"  data,  nor  is  the  lack  of  an  error  flag  a  guarantee  of  "good" 
da !.  .1 . 


3.1.6  RECOMMENDATIONS.  A  positive  neans  of  detect iivj  bad  data  is  to  c(jllect 
three  data  sets  at  the  same  altitode  over  the  same  area,  thus  providioj  a  na^ority 
vote  canparison  to  eliminate  the  ''had*'  data  point.  However,  c.'ns  is  not  expected 
to  be  operationally  practical.  Assiminj  only  one  data  set  is  available,  it  would 
be  appropriate  to  select  an  initial  master  impedance  commensurate  with  ttie 
expected  soil  type.  Should  error  residual  flags  occur,  the  master  i  ii|>edancf 
should  be  iterated  up  or  down  one  soil  type  away  froii  tlie  slave  impedance  limit 
that  was  previously  flagged.  If  fie  residual  TO  error  re.iiains,  't  is  then 
recomi.iended  that  suspect  data  be  renoved  pending  further  data  gathering  aiul 
analysis  of  the  immediate  vicinity  of  data  points  in  question. 
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3.'-^  31;JlCT1V[  2. 


IVo'jf'di;!  Const  ifits. 


3.2.1  DESCRIPTION.  Deteriune  toe  effects  of  chanijes  in  the  vjlaes  of  M.ister 

I  Index  of  Refraction,  and  Vortical  Lapse  Factor  on  the  accuracy  of  the 

vjarp.ije  iiodel  . 

3.2.2  TEST  PROCLDURE.  This  test  was  pe-'Toraed  in  ln^o  parts.  The  first  ,)art 
e<a. lined  the  effects  of  selecting  different  values  for  the  noi.iinal  projra'. 
constants.  These  constants  were  varied  one  at  a  ti.ne  over  the  ratujes  jiven  in 
TaPle  Ii.  For  coi.iparison  purposes  the  theoretical  maxina  and  miniina  are  included. 


TAL„ ;  II 


Constant  Ran<jes  and  No.iinal  Values 


Parai.ietor 

1 

i 

■^est  Range 

F  Nominal 

j”>:a*x""'' 

Maxi.na 

or  l  ipodanco 

O.OOlOSo 

T 

i  0.03 

'  0.036 

■  il.'K) 

‘  1 

0.  )'■; 

>.  Ilf  2  a Tract  ion 

;  1.0302 

1  1.0301 

;  1.330330 

1.3004  ; 

;.  nn4 

ical  l.apse  Factor 

;  0.60 

1  o./o 

'  O.Sb 

0.06  ; 

1.  M 

for  tiio  second  parr,  of  this  o'o,,ective  the  Vortical  Lapse  Factor  was  '"ixed  to 
tilt'  ntjcLiiul  /alue  of  0.3S  at  th-j  inpit  to  the  pr').jraM’,  and  varied  over  t  ie  ranje  ')f 
3.h')  -1.20  in  Overlay  5.  ’’h’.s  t.'valuated  isinj  tiie  tuaeff  ic  ients  at  /aloes  o'" 

/■jrttCil  la.ise  oth-er  thin  t'ne  v.ilae  use  i  to  pr.aduce  ‘he  coe''f ;  c  i  ent  s.  I'l 
add;’  ni,  a  lirief  maiy.is  of  t-ie  effed  of  usinj  coirf'ic  ients  at  vilues  o'"  Index 
of  tract  i O'l  otiu'r  t'nan  v/hat  produced  Mie'.i  was  performed. 

3..'.’.1  ^OGR/W  MODIFICATION.  No  modi  '  i  ;at  i  on  was  nec'tt'ssary  for  part  o  le  of  t  D  s 
0 h.'o'.t , ve.  For  part  two,  tht?  main  pro'jra.a  ivas  anj  le'ited  'oy  the  add. Lion  of 
unit  ,;lo  c  il  1 5  to  Overlay  o  (the  evalaat'.on  overlay).  Each  call  was  aconpanied 
'ey  ,1  'RSi  viliie  of  Vortical  L.jpse  Factor. 

3./. ’.2  DATA  REijUI RLMENTG.  5K  ft  Oasoline  input  data  set. 

3.’/.  5  TEST  RESULTS.  The  statistics  of  tne  first  test  are  consolidated  in  Ta'ale 
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A!5l.'  11 ; 


Mean  Tl)  L-'rors  as  ranclion 
of  [Vjyrain  ConsUinls 


Parameter 

-  *r - 

.■-lean  Error  (usee) 

.  _ 

M’n 

TDA 

Max 

j  N!i  n 

TDB 

Max 

f-Iastor  lai'eddnce 

1 

-.OTIC 

1 

j 

.Ohsh 

1 

-.Oh"/ 

1 

1 

1 

.00'.^1 

Iniiex  of  Re‘'ra;tion 

1 

i 

-.oont 

1 

I 

.001)'' 

) 

^  .0'’'0h 

J 

i 

.Ohio 

Vert’cal  Lapse 

-.OOOt 

1 

.001-' 

'-.ooo: 

1 

.oo;’(' 

For  the  second  les"  'ihe  results  are  given  in  ''dble  Iv'  and  Figures  ?  and  I. 

'  ij  .re  IS  a  ,)1  ot  of  the  .'ean  TD  errors  as  a  ^unction  of  Vertical  Lajise  whil(> 

■  '  'nre  ;  i>  a  ,  1  ot  ef  'he  'can  errjr  in  ♦'ee' .  ’’he  la'(?r  figure  is  a  d'rec!  r’esijl  ! 
)<  •  ee  rst.  !n  'FahFe  IV  the  adj  e  ted  '.lea"  err  )»*  is  the  di‘ference  iie'.an'n  ’he 
eri-ar  j;  each  v.il  i.'  Ve''tic.)l  i  a.'Se  a-',;  ‘he  tsu  err  e-  u'  ,1  '  a;'Sc  'u(.*or  .i* 


’’ABL'  :j 

V,,  ■  re  irs  ac  F".inct'en  e'"  Ve'"  .eal  I  'se 


- » 

'"^4  rrror  [ 

Paw  ‘tcu'i 

^  Ad;js‘'. 

ert ical  Lapse 

1 

4- 

.  _ 

•  rror  (_ft' 

rror 

; 

i 

1 

.40- 

1 

50"  1 

!4  I 

1 

’ '  1 

*  \  7<3 

1  •  *  • 

1 

.  ■’40  ' 

1 

3,M  ' 

1 

;  1.0 

J  ^ 

j  .043 

1 

t 

,  1 

\ 

’75  * 

19 

.  FS 

.000 

« 

f 

.ov  ; 

15., 

1 

0 

.1’ 

'  -.'M3 

J 

-.o/'i  ; 

1 

105  , 

no 

.0'’ 

1  . 

1 

1 

1 

-.23H  1 

3:’'; 

170 

.  1 

j  -..'’14 

j 

i 

-  .  S'l  /  ] 

501 

345 

.20 

1  -.290 

i 

-L 

-.557  ! 

f'f’.h 

533 

I 


3.2.4  EVALUATION  CRITERIA.  Moan  error  is  the  only  evaluation  criteria. 

3.2.5  TEST  EVALUATION.  The  re;iiil  ts  'jiven  in  Table  III  indicate  that  t'ne  choice* 
of  nominal  values  for  tfie  pro'jram  constants  is  not  a  very  critical  decision.  As 
observed  elsev^here  in  this  rejiort,  WAllP  does  not  drive  towards  a  unique  solution. 
The  total  range  of  errors  in  Table  III  is  well  under  10  nanoseconds.  Program 
constants  need  not  represent  the  real  world  values  observed  at  the  time  of  data 
collectici.  Any  differences  appear  as  warpage  and  are  subsequently  corrected  ujr 
in  the  me  lelinj  ()rocess. 

the  s 'cond  test  esseiitially  models  the  ins'-'rtion  in  thi  AN/AR,N-101  of  values 
for  Verti  ;al  Lapse  which  are  different  t'nan  the  one  used  t^  produce  the  v/arpage 
coefficients.  As  can  be  seen  frotn  Table  IV  significant  error  is  introduced,  by 
inconsistent  values  of  the  Vertical  Lapse  Factor.  Since  Vertical  Lapse  is  a 
dynamic  fmction  subject  to  change  both  daily  and  seasonally,  this  test  also 
models  the  effect  of  this  change  on  the  accuracy  of  the  AN/ARN-101.  Although  the 
range  of  Vertical  Lapse  given  in  Table  IV  accurately  models  values  which  may  be 
seen  throughout  ttie  year,  its  effect  on  tlie  AN/ARN-lDl  is  unclear.  This  is 
because  Loran  chains  are  stabilized  rel.itive  to  fixed  monitor  locations.  Changes 
in  the  real  world  parameters  such  as  Master  Impedance,  Index  of  Refraction,  and 
Vortical  i.apse  Factor  are  conpensated  for  within  the  Loran  coverage  area  by  the 
non  i  j)ul  at  ion  of  the  (?mission  delays.  Variations  in  the  program  constants  due  to 
real  world  fluctuation  then  become  transjiarent  to  the  AN/ARN-101.  This  is  true  in 
the  vicinity  of  the  monitor  stations.  It  is  well  beyond  the  scope  of  this  test  to 
attempt  to  determine  the  effects  as  a  function  of  distance  away  from  the  monitor 
1 ocaL ions. 

The  Index  of  Refraction  problem  illustrates  the  problem  of  Chain  variation  and 
control.  Variation  of  the  Index  of  Refraction  throughout  its  reasonable  range  of 
1.0002  to  1.0004  can  introduce  an  error  of  about  100  ft  per  station  at  a  range  of 
1000  Kin.  This  imposes  a  m.)xiinum  ujiper  bound  on  the  position  error  to  about  300 
ft.  Changes  in  the  Index  of  Refraction  cause  a  linear  change  with  distance  in  the 
times  of  arrival  of  the  wave  fronts.  These  clianges  are  compensated  for  at  the 
monitor  locations  by  adjusting  the  onission  delays  by  a  fixed  amount.  This  step 
cliange  is  not  e<pial  ly  valid  everywhere  in  the  coverage  area. 

3.2.b  RECOMMENDATIONS.  The  major  findings  of  this  objective  clearly  state  the 
need  for  consistency  of  variable  definitions  between  WARP  and  the  AN/ARN-lOl. 
Choice  of  variables  is  non-critical  as  long  as  the  values  chosen  allow  the  data  to 
converge  to  acceptable  solutions.  Warpage  coefficients  then  generated  are  valid 
for  real  world  conditions  similar  to  those  that  existed  at  the  time  of  data 
cal  leC  ion. 

It  is  firmly  recommended  that  .idditional  analysis  and  testing  be  undertaken  to 
dc'iermine  tiie  long  term  validity  of  the  warjiage  coefficients  over  seasonal  changes 
as  well  as  during  adver-se  weather  conditions.  This  will  require  a  detailed 
i nvest igat ion  of  tlie  AN/ARN-lOl  coordinate  conversion  process  in  addition  to 
futThei-  restutreh  into  tiie  relevance  of  the  selection  of  constants  to  model 
variables. 
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3.3  OBJECTIVE  3.  Improper  Spheroid  Selection. 

3.3.1  DESCRIPTION.  Determine  the  degradation  of  modeling  accuracy  an  a  function 
of  incorrect  and  inconsistent  spheroid  selection. 

3.3.2  TEST  PROCEDURE.  To  limit  the  scope  of  this  test,  two  cases  were  evaluated 
to  detennine  the" mag ni tude  of  positional  error  resulting  from  selection  of 
inconsistent  or  incorrect  spheroid  (more  precisely  ellipsoid)  earth  nodels.  The 
first  test  was  to  mis. natch  spheroid  coordinate  systems  in  processing  the  data 
within  WARP,  and  the  second  test  was  to  simulate  an  AN/ARN-lOl  operator  error  by 
selecting  an  inconsistent  earth  model  in  the  AN/ARN-lOl. 

Reference  spheroids  are  implicitly  or  explicitly  defined  in  four  separate 
areas  of  the  warpage  correction  process.  "^Iiese  areas  are  t'ne  LORAN  Cnain 
coordinates,  the  LAT-LON  data  pairs,  processes  internal  to  WARP,  and  processes 
within  the  AN/ARN-101.  In  the  AN/ARN-lOl,  any  one  of  seven  splieroids  .nay  be 
selected  and  the  defining  parameters  of  tiiese  earth  .nodels  .ire  given  in  Tal)le  V. 


TABLE  V 

REFERENCE  SPHEROIDS 


REFERENCE 

SPHEROID* 

SEMI -MAJOR  AXIS 

IN  METERS 

INVERSE  OF 

FLAT  TEN  I  NO 

International 

6378388.0000 

29  7. 

Clark  1866 

6378206.4000 

294.978698 

Clark  1880 

6378249.1450 

293.464990 

Everest 

6377276.3452 

300.801699 

Fwjssel 

6377397.1550 

29').  152812 

Austral  1  an 

6378160.0000 

298.25 

■■ICS- 72 

6378135.0110 

298.26 

♦The  reference  source  fo''  these  spheroid  models  is  the  Program  Development 
Speci'^ication  for  Navigational  Computer  of  RF-4C  Digital  Modular  Avionics  System 
A'l/ARN-i:il,  Contract  Number  F19620-76-C-0024,  CB1001-004,  Page  304,  Table  IIP, 
Charting  Spheroids. 


13 


I'l  L'lo  firsL  r.he  Lord.-i  chain  md  warpaijo  data  was  nodelcd  using  the 

survey  '..yst.?u.  The  seven  different  spheroid  iiiodels  were  then  used  in 
je'u'r'ut  ing  warpuge  cot.'f  f  ic  ients.  Positional  errors  were  noted,  thus  simulating 
iiipr  )|)er-  spiieroid  selection. 

tor  test  two,  the  warpage  coefficients  were  generated  using  the  consistent 
sur/  .>y  system  (i.e.,  WGS-72  data  and  spheroid  liiodel ) ,  but  the  seven  different 
spiieruid  uiodels  were  then  used  to  evaluate  positional  error,  thus  si.iulating 
A’lN-ini  operator  selection  of  the  v/rong  earth  nodel . 

3.T. ’.1  ‘PROGRAM  MODIFICATION.  WARP  was  uiodified  to  use  the  Sjiheroid  under  test 
'ns'.eao  of  tiie  fixed  Clark  IBGG  model.  The  only  purpose  of  the  reference  spheroid 
:  ,  in  the  calculation  of  arc  length  distances  froii  tiie  data  points  to  tiie 
s  i'  ions,  ilowevor,  this  calculation  is  done  twice;  first  in  Overlay  3,  used  in 
I  >  generation  of  warpjge  coefficients;  and  second,  in  Overlay  5  used  for 
ev.il  not  iiig  tlio  coefficients. 

!  ir  the  fiesl.  te^it,  WARP  was  ;iodifi,-d  so  tiial  lioth  sptieroids  would 
I  i;  I '..i'loDusly  us“  !.he  values  in  Table  V.  In  the  second  test,  the  first  spheroid 
■/a  1  selected  to  iiu'  tiie  WGS-/3  ’iiodel  whi’e  the  second  was  iterated  tfiroujh  all  the 
iiihle  values. 

3.3.  '.3  OATA  RCQUIREMCNTS.  One  set  of  input  data  and  t'ne  sp'neroid  values  defined 
i  n  I lol 0  V.  ” 

3.3.1  TEST  RESULTS.  The  test  results  are  given  in  Tables  VI  and  VII.  In  Table 
VI  t  ie  spheroids  are  tiie  same  and  in  Ta'ole  VII  only  the  second  siiheroid  is  varied, 
figure  4  is  a  plot  of  the  normal  probability  curves  for  four  different  spheroids 
wii.n  the  Australian,  Clark  1330  and  Cla''k  1366  being  virtually  identical  to  the 
dCS-'T  [)loL. 

3.3.4  EVALUATION  CRITERIA.  Evaluation  criteria  is  the  standard  deviation  and 
hie  in  ori'jr  of  the  total  error  function  given  in  feet  as  shown  in  Tubles  VI  and 


i.J.n  TLSr  EVALiJAriON.  The  results  of  tiie  tests  indicate  that  tiie  accucicy  of 
the  warpuge  correction  is  not  a  function  of  spiieroid  choice  as  long  as  that  choice 
is  onus  istent.  ’’lie  '.vTS-ZC,  Australian,  Clark  1866,  and  Clark  1330,  and  to  a 
lesier  extent  trie  Inn’cnational  appear  to  be  virtually  i  nterciiangeabl  e.  Troii 
I'a'n'ie  VII  ii  is  clear  i.hat  the  Bessel  and  Everest  spheroids  cannot  be  arbitrarily 
-.ills  i  up -j  1  fur  .jiiy  oi'ier  spei’oid.  Ttie  iiajor  hazard  in  the  spheroid  selection 
prov.e-.s  is  to  select  in  the  AN/ARN-lOl  either  of  these  two  Sjhieroids  when  tiie 
war  ,  ijo  coefficients  have  been  generated  using  a  different  one.  Tliere  is  also  a 
P'li ,ii  huzard  if  lie  input  cata  and  the  Loran  chain  are  in  different  rijference 
syste,as.  This  error  soiirce  would  be  very  difficult  to  Pleasure  because  it  would  tie 
a  '.'rang  function  of  position  on  the  earth. 
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TOLL  VI 


I’ositional  Errjr  Due  to  Inconsistent  SoheroiJ  Seli^ction  in  Wnrjui'je 
Coefficient  Cioneration  Process. 


SPHCROID 

MEAN  LRRJR'  (US'EC) 

AVERAGE 

MODEL 

TDA 

TDB 

ERROR  (FT) 

INTERNATIONAL 

-.DPOC 

n'M 

133.20 

CLAtL  ISRG 

-.0302 

.0017 

133.44 

LVERLSr 

.  "000 

134.49 

3lssi:l 

.'T'Rl 

."0C3 

:  13.34 

AUSTRALIAN 

.0002 

."018 

133.44 

NGS  !?- 

.0003 

.0  018 

133.48 

CLARK  1S66 

.0001 

mi  • 

133. 4‘) 

L. 

TABL!  VI  I 

. 

pQsi tionai 

r_r»'jr 

Due  To  Incor'-ect  3n'neroiil  So’ec 

■  ion  in  the  AM/A'Li-l  'l . 

SPHEROl  I 

r 

MEAN  L  ;R0R'  C'JS.'C) 

'■  AVERAGE 

MODEL 

1 

TDA 

TDB 

ERROR  (FT) 

INTLRNA'^IONAL 

- 

1 

i 

.0003 

.  0018 

133.48 

CLARK  13  LU 

-.0030 

_  noi,! 

133. An 

CVTRTST 

. 

( 

.l'^04 

.^313 

296.4b 

CLSSLl. 

.I4bn 

.  158L 

327.02 

AJSmIALIAN 

1 

o 

-.  0007 

134.J3b 

NGS  /:’ 

1 

1 

-.0231 

.002? 

134.0b 

C..ARK  13n() 

1 

1 

1 

L- 

-.044', 

-. 0/01 

L _ _  , 

153./') 

1  . . 

WGS-72 

INTERNATIONAL 


Figure  Spheroid  Probability  Functions. 


3.3.6  RECOMMENDATIONS.  WARP  should  he  'HodirieJ  to  perfurn  the  distnco 
odl  C'jl  at  i  ons  only  once.  This  will  >"equire  t'nat  the  distances  he  saved  in  cjiaon 
storage  for  use  by  additional  overlays.  In  addition,  the  constants  l■.■.-pre^,c'nt  i  ng 
all  the  various  spheroids  should  be  added  to  the  jirograin.  The  jirogrj  i  op-oMtoi' 
should  be  able  to  select  any  of  the  available  spheroids,  '^hese  changes  w’ M  s.iate 
WARP  consistent  with  AN/ARN-101  operation. 

Additional  analg/sis  needs  to  lie  done  on  tiie  re  gui  reneni.  for  the  A'l/A!tN-I '11  ti 
nodel  seven  different  Sjihoroids.  I'able  VI  impl  ios  that  tin'  choice*  of  ■,  liu*--:) i.l  ;  , 
irrelevant.  the  reference  spiit'roid  is  only  'jsed  '  n  i:o!i]'ute  ar-c  L'-ipii',  v'en 

t'ie  data  points  and  trie  i'.'lRAN  st.i'  ions,  trien  seleotalil  c’  spiirer-.o  i>iv  .'ojld  In- 
el  iiii  noted  in  favor  of  one  sptieroid.  This  woulri  p'-jbably  increase  the  -ipp-ireo' 
warpage  but  it  would  be  a  sel f-correct  i  ng  process.  As  long  as  WAR”  aril  t hr 
A'l/ARN-lOl  are  consistont  there  should  be  no  degradat  ion  in  the  overi;’  no  ie’ 
accuracy.  The  operation  of  the  AN/ARN-IOl  and  WAR'’  would  ne  s  impl  *5 i-ad. 
Additionally,  some  savings  in  AN/AkN-lOl  prograui  size  mi  iixecition  ■  le  "n.’t!  be 
real i zed. 

The  Defense  Mapping  Agency  Aerospace  Center  (T-'AAC)  has  the  initial 
responsibil  i  ty  to  exorcise  WARP  and  gener.-te  .varpago  Ciieff  icients.  D'-'AAC  has  thi* 
capability  to  transl.ite  LAT/L9Ns  '‘ron  one  coordinate  system  to  anot'nm  i  is 
not  particularly  important  to  fix  the  cho’C'a  of  referonce  'oefor'e  t'ne  inpd.  d.it  i  ' 
gatriered.  It  is,  howevo*',  most  important  I'nat  the  user  speci'^y  the  sp'ncroidmi) 
that  is  desired  fo>'  iiis  operational  area.  For  tiie  best  -m.cirricy  the  A'!/A'-::;-l ')l 
'■lust  use  the  sane  spheroid  that  w-is  used  •'o  generate  the  coef i c  ion'i  *.. 
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3.4  OBJECTIVE  4. 


!n,;ut  Daia  Der.sity/Si/e  of  PCA. 


3.4.1  DESCRIPTION.  Investigate  the  degi-adat i on  of  )josition  accuracy  as  a 
function  of  !n;jut  data  density. 

3.4.;?  TEST  PROCEDURE.  Input  data  sets  were  generated  at  the  following  densities: 
5nin  sg,  I'rtnni  x  5n,a,  lOraii  sq,  IBnin  x  Snin,  and  ISniii  sq.  The  data  set  at  the  Sn.n 
density  was  the  evaluator  data  set.  The  other  four  data  sets  were 
input  to  WARP  and  warpcige  coefficients  generated.  Eacli  set  of  coefficients  was 
retained  and  evaluated  using  the  evaluation  data  set.  This  consisted  of  rerjtining 
prugraia  Overlays  1,  2,  3,  and  5  with  the  evaluator  data  set  as  input  to  the  basic 
prugra.a.  This  enabled  evaluation  of  coefficients  at  a  data  point  density  greater 
than  the  input  point  density. 

3.4.2. 1  PROGRAM  IIODIFICATIQN.  Two  progra.i  iuoji  f  ications  were  uade.  The  first 
c'nange  .vas  to  add  ihe  abil  i  ty  to  attach  two  input  files,  with  different  nd'ses,  to 
the  natch  Job  control  file.  The  second  oodif ication  was  made  to  '4ARP  itself. 

Spec i i cal  ly ,  Overlay  0,  liie  pro-jram  executive  routine,  was  extended  after  the 
original  final  overlay  call  to  redefine  the  inpot  data  file  to  the  second  data 
file  iti.ached  by  the  Batch  Job  Control  f:le.  Then  all  overlays  except  Overlay  4 
were  called  again.  Since  the  coefficients  v/e''e  generated  by  Overlay  4  and 
reiiuii'ied  in  coiwnon  storage,  this  had  the  effect  of  evaluating  the  second  data  set 
isinj  the  coefficients  as  generated  by  p’'Ocessing  the  first  data  set. 

3.4. 2. 2  DATA  REQUIREMENTS.  Eron  a  mast?r  data  set,  an  eval uator  data  set  of 
points  at  a  density  of  approxima .ely  1  point  per  onn  square  v/as  generated.  The 
ind. .‘pendent  data  sets  listed  above  were  also  created.  Less  than  10%  common  points 
were'  allowed  botwO'-'n  the  test  and  evaluation  data  sets. 

1.4.1  JL^Sr  ^RESULTS.  The  test  r..'S.:lts  are  given  in  Table  VIII.  The  asterisks  in 

(.he  ible  indicate  errors  greater  than  allowed  by  the  cOnt-witer  progra>ii. 

I. 4.4  LVAL'JATION  CRITERIA.  The  ev.iluation  criteria  is  t!ie  change  in  the 
coiiipojite  nean”errdr  Tunction  in  feet  for  the  test  data  sets  as  conpared  t.)  the 
error  of  (.ne  oval  ntor  data  set. 

3.4.2  TaST  EVAlUATIOM.  Frou  Taale  VIII,  it  is  concluded  shat  data  gathered  at  a 
lenu.  ,y  of  less  f.h.an  1  p.aint  eve-'y  5n..'.  dees  not  adequately  model  the  coverage 
ar./i.  ’"here  is  also  so  ;e  oupirical  data  which  indicates  that  'With  a  larger  pri.ae 
■  If.'..  !  les'S.:r  d.jta  density  nay  be  tolera'am?.  Unfortunately ,  no  data  over  such  a 
lurji;  aro.i  exists  but  so  ae  observ.jtions  .aay  be  'aade. 

2  vitour  ,;lo,s  of  tiru  various  data  sei.s  show  that  'jniform  data  cover.age  is  very 
;  ;  '  ir'tun:..  Void  jfe.i-,  Jo  not  constrain  'he  ■aode'ied  s urface  and  many  may  ictually 
,r  )i;,)  ..’  Vi?ry  radical  oehavio’'  of  the  uiodel .  Evaluation  of  the  node!  only  at  tiie 

J. ,.’  I  paints  wi  1  ’  -lOt  dn.ect  this. 


TABLt  V!  ;  1 
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The  data  at  Eglin  AFB  was  gathered  over  a  region  about  40  X  bOnn;.  At  t ne 
greterred  density  of  1  point/Snni  square,  tin's  gave  AO  points  organi pi;  1  i  <  ’.I,  A 
prime  area  of  lOOnm  per  side  will  have  40C  points  at  t'ne  same  density  urganipe.i  ?.() 
K  20.  Doubling  tFie  interval  between  the  data  points  will  give  4  x  A  and  lil  X  li) 
points  ove’"  the  respective  coverage  areas.  Although  the  two  areas  di'^far  greatly 
in  size,  this  distinction  is  nieaningless  within  WARP  and  within  the  AN/ARN-lOi 
system  as  well.  Every  ori/ae  area  is  noraal  i zeu'  into  a  region  of  fixed  diiiens'ons 
regardless  o^  its  s'ze.  Since  the  vwrpage  moiie'  is  preordained  to  be  of  order 
four  it  would  see;)  t'nat  mathematical  ly  the  numduer  of  data  points  constraining  IFie 
tv/o  dimensions  would  be  vc'y  important.  '-;ov;eve’',  furt'ier  Mivestigation  of  the 
nathemat  ical  requ  i  renents  verses  t'ne  geop'-'ysical  realities  of  t'ne  ear'.h  is 
considered  to  'ue  beyond  ‘.'ne  sco;)e  of  this  test  effort. 

3.4.6  RECOMMENDATIONS.  Although  there  is  so;iie  ov'dence  to  the  contrary,  the  lac'x 
of  data  over  a  large  area  forces  the  conclusion  that  data  should  be  gathered  on 
Slim  centers.  Since  void  areas  ca  ised  miioun'Ied  errors  to  occur,  it  is  reco'iiuondi'il 
t'nat  WARP  be  nodified  to  fill  in  empty  cells  wi  t'n  roj^onahle  i  npedanrr  valijs. 
fhis  will  tend  to  constrain  the  wirpage  polyno.,'ial  to  reasonable  warpage 
correction  values  and  ‘^  /■"tber  a^d  in  "smo  ;t'ii nj"  tiie  nodeling  process. 


OBJECTIVE  5.  T.D.  Prediction  in  Void  Areas. 

3.S.I  OESCRIPTI ON.  Determine  the  ahility  of  the  proyram  to  predict  time 
diTf. 'fences  in  regions  of  the  prime  are.-,  that  are  devoid  or  nearly  devoid  of  inpet 
points. 

-  TEST  PROCEDURE.  N/A 
B.B.'M  PROGRAM  MODIFICATION.  N/A 
3.J.P.?  DATA  REQUIREMENTS.  N/A 
3.3.3  TEST  RESULTS.  N/A 
3 . •■) . 4  EVALUATION  CRITERIA.  N/ A 

3.3.‘i  TEST  EVALUATION.  At  the  time  the  method  of  test  was  written  it  was  felt 
tnal,  .1  prediction  ca.nJbility  might  exist  within  the  program.  The  results  of 
Objective  4  and  an  analysis  of  the  comp  iter  program  clearly  shows  that  no  sndi 
c ipab i 1 i ty  exists. 

3.0.n  RECOMMENDATIONS.  See  Objective  4. 


3.6  OBJECTIVE  6.  Altitude  Time  'difference  Effects. 


3.6.1  DESCRIPTION.  Determine  i .J  data  set  jat'nered  at  one  altitude  can  jiMier.iti 
coefficients  which  are  usaiile  at  other  alt  tudes. 

3.6. ?  TEST  PROCEDURE.  Four  sets  of  data  were  gatiiered  over  the  Eglin  '■’rime 
Coverage  Area  at  the  following  approximate  altitudes:  1000,  6000,  10000,  and 
16000  ft  MSL.  These  sets  were  independently  edited  to  a  density  i)f  one  point  per 
five  nautical  miles  sguare.  The  evaluation  of  the  liata  pron'eeded  ui.iilirly  tj  I  lu 
■iiet'nod  used  for  the  Input  Point  Density,  Objective  T.  That,  is,  each  il.Da 
generated  a  warpage  model  vvhich  is  completely  descrilaed  mathematical  ly  by  tiie 
warpage  coefficients  and  ti)e  area  boundaries.  All  ^our  data  sets  v/oro  ju-jcessecl 
against  each  of  trie  four  nodels  and  the  ofects  on  eac'i  was  i.ieas  ired.  These 
sixteen  combinations  of  data  and  warpage  models  gave  a  co.njirehensi ve  view  of  the 
alt  i  tide  effects  for  the  prime  coverage  ar-aa  under  test. 

3.6.  ?.l  PROGRAM  MOOIFICAf ION.  This  test  objective  used  t'ne  same  version  of  ’JAPI’ 
as  did  Objective  A. 

3.6. ?.?  DATA  REQUIREMENTS.  Pour  sets  of  '.AT-LDN/TD  data  jiairs  yathere!  within 
t'ne  following  altitude  boundaries:  0-2500  ft,  ?600-750'J  ft,  7500-l?6'lh  ft,  .md 
I2600-17500  ft. 

TEST  RESULTS,  fable  IX  shows  the  average  positional  error  and  riurhial 
standard  deviation  of  the  data  as  well  as  the  ()narf.ile  boundaries,  tai)1e  is 

nrgarpzed  in  four  parts.  !n  eacii  part  a  different  coef f  ic  ienl-genera' i  Hj  input 
data  set  is  held  constant  and  all  '^oir  data  sets  .we  evaluated  against  't. 

Table  IX  does  not  include  data  collection  error  as  part  of  trie  dat.'. 


TABLE  IX 


'.IRAN  WARPAGE  POSITIONAL  ERROR 


a.  WARPAGE  MODEL  GENERATED  FROM  l'<  INPUT  DATA 
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Ii  I,  -..‘(in  '.'idt  ■  n-'  si.iooLti  i  ri'j  fjnc!:ion  of  the  noilt.'l  i  o  j  process  co.iipl  er.ents  the 
s  ioo'h.nj  function  of  nov.nj  jway  ^rj:i  'he  earth's  surface.  It  is  probably  not 

to  begin  v^i  t.h  saootiier  high  le/el  data  and  try  to  roughen  it  sufficiently 
"  )  iescribe  warpage  at  lower  altitudes  uecause  there  .light  be  local  ano.uilies 
wh'ch  are  evident  at  lower  altitudes  only.  Metal  structures,  jiower  lines,  and 
fcnc.'-,  producH  effects,  for  exa..)j)1e,  that  may  be  important  at  IK  feet  but  nrt  hK 


l.n.n  ■If COMMlNDAflONS.  In  flat  terrain,  one  or  two  low  altitude  (bita  sets  should 
t  ..r  ely  .describe  the  warpaje  as  a  function  of  altitude.  Additional  levels  of 
li'  I  tend  to  be  f  unci  >  ona  1  ly  roJunilant  lut  can  be  very  valuable  for  data 
v.i!;  I,c  ion  ijuruouc's.  ’'hree  levels  of  d.ita  are  sufficient  to  provide  for  a 

ir  I vote  .10  ms  .if  y.ilidation.  Validation  is  very  important  because.*  of  the 
.Mt''.‘..io  d.t'  ic.ilty  of  sepe.ratinj  .iieasur  Mont  errjr  fro.ii  warpage.  Data  gathered  at 
I  i.'/’n  ultitides  IS  ..lore  descriptive  of  the  total  vertical  profile  than  data 
ji''i.'r;d  ,|L  higher  altitudes. 

'o  in^.ire  '.b.‘  .lost  pr._*cise  .nav'gati.in  solutions,  data  shauld  be  gathera.l  ,3t 
i‘  f'o.r  jl  t  i  .  i;1e  levels  specified.  Realizin.j  thai,  .ivailablc  resources  .i.iy 
;e;'  c.e  that  fewer  nltit.ides  bo  flown,  it  is  strongly  recomiiiendod  that  data  be 
c  d’..;cted  at  tne  s.ime  iltitide  at  whicli  s.jliseguent  operational  nissirns  are  to  be 


3.7  OBJECTIVE  7.  ;  . 

3.7.1  DESCRIPTION.  ,  '.'i-V  Ounddry 

condi t 1 ons. 

3.7.2  TEST  PROCEDURE.  '  'w  sf ,.oint5 

consisting  of  i)d'''o.'  •  .  .  .  ,  ,  r  •  d'stdnce  .icross  the 

prime  area  I)oijndd’'y.  '' >j'  >  '  •  •  ' '  ■  ,,oint  tds  fron  tiie 

prime  area  coef  ic  ler' ■,  .m'  •  •  .  ••  j  •  mp  cdlcilatr-l  constant 

impedances  of  the  appro^-’f  •  ^ v  ■;  .  ■■^‘■erence  i)et.n3L‘n  tlie 

interior  and  exterior  ’'Ds  ^a  -  ‘  •  n.’  rplj-.i/e  valjo  of  tiie  position 

jump  at  the  boundary. 

3. 7. 2.1  PROGRAM  MODIFICATION.  ')ve’'ldy  S  W"'-’*’  ,vds  al  tered  to  create  ten  pairs 
of  data  points  equally  space-.^  aionj  eacn  o''  'he  '"om  prime  area  boundaries.  Each 
data  pair  consisted  of  a  point  m  ‘.he  in‘.e’'ior  of  the  prime  area  within  300  ft  of 
trie  boundary  and  a  correspond' ng  point  the  same  1' stance  across  the  boundary  into 
a  secondary  area.  The  approx-. laate  relative  locat’ons  of  the  data  points  are 
illustrated  in  Figure  13.  Additiona’  progra-.i  lodh  "i  cat  ions  were  implemented  to 
compute  TDs  as  follows: 

TDA  =  TDA  (exterior)  -  TDA  (interior) 

T06  =  TDB  (exterior)  -  TOB  (interior) 

3. 7. 2. 2  DATA  REQUIREMENTS.  Low  altitude  data  set  processed  by  normal  Overlays 
1-4  and  a  "modified  Overlay  5. 

3.7.3  TEST  RESULTS.  Test  data  is  compiled  in  Tables  K  through  XI 1 1.  Each  table 
is  organized  by  secondary  area  number  and  contains  TDs,  northerly  and  easterly 
position  jumps  in  feet,  and  absolute  magnitude  of  the  relative  position  jumps. 
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Figure  13.  Boundary  Jump  Errors. 


TABLE  XII 


AREA  4 

SQUNOARV 

CONDITIONS 

IDA 

TOB 

N  JUMP 

'  E  JUMP 

MAG 

PAIR 

(usee) 

(usee) 

.  (ft)  . 

(ft)  . 

(ft) 

4.1 

.172 

.332 

543 

68 

547 

4.2 

.060 

.268 

444 

16 

444 

4.3 

-.078 

.189 

328 

-  48 

332 

4.4 

-.232 

.108 

214 

-121 

245 

4.5 

-.366 

.041 

126 

-189 

227 

4.6 

1.277 

.172 

75 

625 

629 

4.7 

1.456 

.265 

182 

703 

727 

4. a 

1.422 

.226 

145 

694 

709 

4.'J 

1.270 

.134 

70 

636 

640 

4.10 

1.081 

.050 

13 

560 

560 

TABLE 

XIII 

TOnT' 

AREA  5 

BOUNDARY 

CONDITIONS 

fO'A" 

fOB 

N  JUMP 

’  E  JUMP 

MAG 

PAIR 

(usee) 

(usee) 

(ft) 

(ft) 

(ft) 

5.1 

-  .595 

.418 

641 

-353 

732 

5.2 

-  .821 

.060 

163 

-419 

450 

5.3 

-1.027 

-.163 

-150 

-497 

519 

5.4 

-1.194 

-.220 

-225 

-575 

618 

5.5 

-1.315 

-.098 

-  4 

-647 

647 

5.6 

.002 

.412 

747 

-  26 

747 

5.7 

.016 

.378 

604 

-  13 

604 

5.3 

.015 

.300 

530 

-  11 

530 

5.9 

-  .010 

.274 

482 

-  21 

483 

5.10 

-  .065 

.243 

431 

-  46 

433 
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3.7.4  EVALUATION  CRITERIA.  The  key  parai.'eter  evaluated  was  the  'tiagnitude 
position  error  shown  in  tables  X  through  XIII.  This  position  error  is  the  result 
of  the  northerly  and  easterly  errors  and  i^epresents  the  instantaneous  change  in 
Loran  position  observed  by  the  AN/AP.N-101  as  the  prime  area  bc.'ndary  limit  is 
crossed. 


3.7.5  TEST  EVALUATION.  Although  the  test  results  shown  in  Tables  X  through  XIII 
are  val id  Tor  the  data  gathered  at  Eglin  AFB,  they  are  generally  applicable  to 
other  operational  areas  as  well.  At  Eglin  there  is  a  wide  variation  in  api’arent 
impedance  ranging  froiT,  near  saltwater  along  the  southern  boundary  to  poor  soM  in 
the  north.  This  larga  variation  in  impedance  over  the  region  o"  interest  s 
probably  typical  of  virtually  all  areas  of  operation  for  the  AN/ARN-IOI,  with  Hie 
possible  exception  of  rugged  coastal  mountainous  terrain. 

As  indicated  by  the  data  in  Objective  3,  the  uncorrected  warpage  contributes 
about  1500  ft  position  error.  The  price  of  enjoying  a  200  ft  Loran  navigation 
accuracy  in  the  prime  area  must  be  paid  at  the  boundary  crossing  because  of  the 
switch  from  a  precise  to  an  average  impedance.  Therefore,  it  is  reasonable  to 
predict  worst  case  position  jumps  on  the  order  of  1000  ft  or  more  for  almost  all 
prime  area  boundaries.  This  jump  is  most  severely  reflected  in  the  LORAN/LORAN 
mode  of  the  AN/ARN-101  operation.  The  Kalman  filtering  action  in  the  other 
AN/ARN-101  modes  will  not  allow  i nstantaneous  position  jumps. 

Because  WARP  is  designed  to  be  used  with  as  many  as  four  different  data  sets, 
four  slightly  different  sets  of  constants  and  coefficients  are  probable.  Each  set 
would  result  in  somewhat  different  boundary  conditions.  Although  the  AN/ARN-101 
can  handle  four  altitude  sets  of  polynomial  coefficients,  only  one  constant 
impedance  is  stored  for  each  secondary  area.  The  choice  of  secondary  impedance 
values  is  relatively  unimportant  since  the  altitude  effects  of  warpage  are  at 
least  an  order  of  magnitude  less  than  the  boundary  effect  (as  shown  by  Eglin 
data).  For  th. s  reason,  any  secoridary  area  constant  impedance  generate^!  by  a 
valid  set  of  data  will  suffice. 

Although  not  directly  related  to  the  quantification  of  the  boundary 
conditions,  several  program  limitations  were  noted.  At  least  one  real  data  point 
must  reside  in  any  cell  adjacent  to  a  secondary  area  border  in  order  for  the 
program  to  calculate  an  impedance  for  the  area.  Also,  data  points  in  a  secondary 
area  do  not  contribute  to  the  calculated  secondary  area  impedance.  All  data 
points  are  used  in  the  calculation  of  the  PCA  polyno'nial  coefficients  whether  they 
are  interior  or  exterior  to  the  PCA. 

3.7.6  RECOMMENDATIONS.  Several  imp'*ovements  to  reduce  the  boundary  effect  can  be 
iiade.  One  is  to  change  WARP  and  the  other  is  to  carefully  select  the  secondary 
area  constant  impedances  to  be  inserted  in  the  DMAS  computer. 
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As  currently  imp! anenteci ,  WARP  searches  the  perimeter  cells  of  the  prime  area. 
An  average  impedance  for  each  secondary  area  is  calculated  from  the  data  points 
resident  in  the  perimeter  cells.  However,  the  program  could  be  modified  to 
include  data  points  in  the  secondary  areas  near  the  border  (i.e.,  5nm  or  less). 

If  the  additional  points  are  well  distributed,  a  better  average  impedance  should 
result.  No  special  effort  to  gather  data  outside  a  prime  area  should  be  made  for 
titis  purpose  as  the  improvement  would  be  small.  However,  prime  areas  are  not 
expected  to  be  isolated  and  if  data  earmarked  for  one  prime  area  is  available  for 
use  in  a  different  PCA  secondary  zone  it  could  be  used.  This  approach  is  flawed 
by  requiring  engineering  judgement  and  data  base  management  not  currently 
envisioned  for  WARP  operation.  WARP  should  be  modified  to  detect,  reject  and  list 
all  data  points  outside  the  PCA.  This  would  simplify  WARP  and  provide  an 
additional  error  check  on  the  input. 


There  does  exist  a  cumbersome  but  very  effective  method  of  suppressing  the 
boundary  juiiip.  It  requires,  however,  A  PRIORI  knowledge  of  where  the  boundary 
crossing  will  occur.  If  this  is  known  then  the  impedance  of  the  data  point 
nearest  the  crossing  can  be  substituted  for  the  calculated  secondary  impedance. 

The  necessary  data  is  available  as  part  of  the  current  program  output  but  it  isn't 
foraatted  to  this  end.  It  is  recommended  that  the  WARP  program  output  be  modified 
to  add  a  table  or  plot  of  the  boundary  impedances  so  that  the  users  have  the 
freedoin  to  select  the  values  that  best  suit  their  needs.  How  such  a  selection 
process  would  be  implanented  in  the  field  is  unknown,  but  the  proposed  AN/ARN-101 
Mission  Data  Transfer  System  is  powerful  enough  that  this  capability  could  indeed 
be  added. 
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3.8  OBJECTIVE  8.  Average  and  Saltwater  Impedances. 

3.8.1  DESCRIPTION.  Determine  and  canpare  the  positional  accuracy  within  the 
prime  area  using  saltvvater  and  average  impedances  instead  of  the  polynoiiiial 
warpage  model  impedances. 

3.8.2  TEST  PROCEDURE.  Three  sets  of  constant  impedances  were  evaluated  to 

determine  their  relative  merit  in  correcting  for  Loran  warpage  within  the  prime 
area.  This  was  accomplished  by  modifying  WARP  to  use  selected  values  of 
impedances  in  lieu  of  the  warpage  model  impedances  calculated  from  the  input  data 
set.  Table  XIV  compares  the  various  impedance  values  used  for  the  naster  and  two 
slaves  (Z|T,,  Z^,  ). 


TABLE  XIV 


CONSTANT  IMPEDANCE  VALUES 


nMe 

^m 

Z2 

Saltwater 

.001055 

.001055 

.001055 

Mil  1 i ngton 

.040 

.028500 

.014410 

Best  Ave  D 

IK 

.040 

.0431456052 

.0233970566 

Best  Ave  P 

5K 

.040 

.0439589757 

.0233413805 

Best  Ave  !? 

lOK 

.040 

.0450740286 

.0230334424 

Best  Ave  (? 

15K 

.040 

.0431774379 

.0223470084 

The  Saltwater  impedance  model  represents  an  ideal  case  which  assumes  the 
entire  Loran  coverage  area  is  ho:nogeneous  with  an  equivalent  impedance  to  that  of 
sal twater. 

Millington's  impedances  were  calculated  by  a  manual  method  of  fractional 
parts.  Basically,  each  arc  length  distance  from  the  data  point  to  the  stations 
was  subdivided  into  elements  of  homogeneous  impedance.  After  each  fractional 
length  was  assigned  an  impedance  value,  they  were  weighted  by  their  fractional 
contribution  to  the  total  propagation  path  under  consideration.  The  weighted 
impedances  were  simply  totaled  to  obtain  the  final  impedance  value.  The  data 
point  used  for  these  calculations  was  the  center  of  the  prime  area. 

The  four  sets  of  Best  Ave  impedances  were  computed  by  WARP  using  the  baseline 
input  data  sets  for  each  altitude.  These  values  are  identical  to  the  average  of 
the  impedances  represented  by  the  coefficients. 
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PROGRAM  MODIFICATION.  Several  .nodi  ^ications  to  Overlay  5  were  .idOc.  "oe 
first  change  caused  WARP  to  canpote  the  average  impedance  for  eacii  slave  station 
v;illiin  the  jirime  area  region.  The  second  modification  forced  Subroutine  Wave  to 
use  selected  impedances  instead  of  calculating  them  from  the  coefficients. 

3.8. P.P  DATA  REQUIREMENTS.  Four  sets  of  altitude  dependent  input  data  and 
impedance  values  defined  by  Table  XIV. 

3.8.3  ^FST  RESULTS.  Table  XV  shows  the  results  obtained  by  using  the  three  types 
of  constant  impedances.  Fach  altitude  set  is  compared  to  the  position  errors 
predicted  by  the  coefficient  model,  figure  11  is  a  plot  of  the  nornal  probability 
functions  using  the  IK  ft  altitude  portion  of  the  data  in  Table  XV. 

3.8.4  CVALUATION  CRI'IlRIA.  A1  '  input  data  points  were  required  to  iiave  real 
world  values  of  cal cul ated  impedance  (.001055  Z  .08).  For  comparison 
purposes  only,  the  WARP  model  generated  at  each  altitude  was  used  as  the 

s  tandard. 

3.8.5  TEST  EVALUATION.  Of  the  three  sets  of  constant  ii.ipedances  evaluated,  the 
Rest  Ave  set  provided  the  best  vvarpage  .iiodeling.  Unfortunately,  this  requires 
t'nat  correlated  TD/LAF-^ON  data  be  gathered  and  processed  to  obtain  a  set  of 
impedances  over  a  geographical  area  of  interest.  This  data  set  does,  however, 
indicate  a  practical  limit  of  tne  accuracy  of  the  Millington  method  in  estimating 
constant  i'npedances. 

As  can  !)e  siien  fron  Table  XV,  no  set  of  constant  impedances  corrects  for 
warp, age  with  any  degree  of  precision.  Only  an  approximate  50%  improvement  is 
noted  fruai  tiie  s.iltwater  model  to  the  best  average.  The  large  mean  error  is 
attributed  to  i ndeiiendentl y  averaging  the  slave  impedances  instead  of  ascertaining 
values  which  aiiniinize  the  total  error. 

Mii’ington's  estimation  procedure  will  not  greatly  improve  warpage  correction 
as  coiiipared  to  the  saltwater  impedance.  In  terms  of  the  AN/ARN-101  operation,  the 
a!)ility  to  select  constant  impedances  for  the  prime  area  in  lieu  of  saltwater  or 
WARP  model  generated  impedances  is  virtually  redundant. 

3.3.5  RECOMMENDATIONS.  Should  it  be  necessary  to  use  constant  impedances  for  the 
priine  aroc;,  two  methods  have  practical  utility.  The  first  .nethod  is  to  acquire 
dciL.i  .ising  a  mobile  or  portable  receiver  in  the  field.  Impedances  can  be  hand 
c.ilculaled  .ind  averaged  using  well  defined  techniques.  The  most  difficult 
computat ional  load  is  obtaining  or  calculating  the  arc  length  distance  from  the 
data  point  to  the  Loran  stations- 


■Another  methiul  of  usiny  cansfont  i'n|)e'Jci'tc<‘s  's  ‘  )  "ooiti*  f'*,"  t  ne  i.i’n. 

o'  0  known  locotion.  ’'o  ho  this  *he  AN/ARN-101  •'eceive''  ■u.st  o*'  i'  'hr  Ir-,'--" 
s,'Ot  with  well  Nnjwn  coorhinates  (i.e.,  second  ('•der  s  irwey  ,ioi!i'. '*  ofis’ 

’  ■);)eddnces  are  inserted  into  the  <AN/ARN-1 ,  the  como  iter  will  jse  'Jie  i 
coordinate  conversion  ;'rucess  to  'jenerate  a  LA"^-'  'N.  !'>y  -teratMij  i  he 
slave  impedances,  the  coneuted  position  car  l)e  ladi  '■■  c-ewerje  to  the  i  n  iw” 
location.  No  A  PRIORI  toran  know'edje  is  "ecessary  and  the  iMpedanco's  jm'- e 
will  he  valid  in  ttie  vicinity  of  the  data  point.  '^he  aajor  hacard  to  ‘ '  s 
approach  is  possible  warpaije  doe  ‘o  strictly  local  pnenoaena.  "^ins  type  of 
warpage  may  be  caused  by  netal  buildings,  fences,  telejihone  lines,  iiov/er  '  i>-, 
e'c.  a  point  fit  is  done  in  an  area  clear  of  locally  disturbing  infl  i.‘n,  , 

usable  results  should  be  obtained. 
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‘^BJlJTIVL  9.  PnoiiJ  Rocunna  . 'janci.*  Ddtd  Validation. 

I.  *.',  'll SCRIPTION.  Validate  the  real  world  data  collection  and  Loran  wur,;aije 
Cix'i'l'.,:  ion  process  n-i'tij  p^ioto  reconna  i  sance  data. 

3.'1.  ’  PROCtPdRL.  Js'.ivj  liie  OK  baseline  data  as  a  reference,  co-npare  the 

tutv  ■  j  j '  ional  acciirKy  wi'nin  the  f^jlin  ^CA  using  data  retrieved  fro;ii  proposed 
h'-dl  pnoto  r-e:  innj ;  5,  j  ,.,e  '’I:  gut  missions.  Compare  reconna  i  ssance  data  with  tiwie 
.lo’-ja r- ;da r  ' phoiotneodol  ;  te  and  RT-AC  ;  nstrumentation  data. 

l.d.M  ’RCGRAr  ‘-lODinCATIdr;.  None. 


i.  '  ^A"A  Rl'C -IRl  ^hcto  reconna i stance  data  and  time  aierged 

r  1  .  oi;  ,)t'\eod.;';l  1  itvi  R-'-AC  1  ns t ru  :ie'd  a t  on  data. 


l.  i.-  _7  RLSULTb.  Alt'iowg’h  '^ijU  Jissions  at  Ljlin  were  flown  to  caaplete 

■ ■>  ;•  v.':,  tno  reconr.a '  s  i  r  ce  film  was  evaluated  by  DMAAC  and  deemodi  anusable 

'  ;  ,  r,M  ic  ;  .  on. 

3.'\A  n.ALdA^’ON  CRITLRIA.  None. 

3 . ') .  h  T:.,  S"  -VAu  None . 

l.^^n  3..:,VtMaNbA:;0NS.  None. 


i.r!  OBJECTIVE  10.  Software  Program  Documentat i on. 

3.10.1  OESCRIPTION.  Identify  computer  software  deficiencies  and  program 
modifications  to  provide  a  .nore  user-orient(?d  operation. 

3.10.3  TEST  PROCEOURE.  The  ident i f icat ion  of  program  deficiencies  and  necoss. 
user  docu'nentation  resulted  in  a  Statement  of  Work  given  as  Appendix  B. 

3.10.3.1  PROGRAM  MODIFICATION.  As  defined  in  Appendix  S. 

3.10.2.2  DATA  REQUIREMENTS.  Not  applicable. 

3.10.3  TEST  RESULTS.  The  tasks  identified  in  Appendix  B  were  accO'nijl  i  shed  liy 
An-nament  Systems,  Inc.  through  a  contract  effort  by  the  Naval  Weapons  Center, 
China  Lake,  California. 

the  results  of  this  Gffo''t  include  a  revised  source  program,  a  User's  Manj 
and  an  Analyst's  Manual. 

3.10.4  EVALUATION  CRITERIA.  Not  applicable. 

3.10.5  TEST  EVALUATION.  Not  applicable. 

3.10.6  RECOMMENDATIONS.  Final  resolution  is  required  to  estalilish  tiw 
resjionsible  organization  for  software  operation  and  maintenance  of  this  prjgra 


4.0  COriCLUS  I  jNS.  Thij  rjoioiidl^  in  testing  the  '.■iar.jijo  Coe  ■ ;  c  i  ent  Cen^.-rat  ’  on 
nroyrj.a  has  been  to  o^jecify  the  program  1  imi  tat’ons,  identify  poss^ole  p'itfalls  or 
prob’ioa  areas,  and  to  gauge  the  level  oi  expertise  necessary  to  provide 
operational  grid  war,;age  support.  The  major  conclusions  of  this  test  d':Livity 
utilizing  Cglin  CCA  data  are  briefly  suinmari ze*!  in  this  section. 

1.  Editing  inpiit  ddta  and  selecting  an  dppr.jpriate  iiaster  stotiiin 
iripeJance  requires  a  high  level  of  program  interaction  and  demands  a  high 
degree  of  Loran  expiertise  (i.e.  program  flags  do  not  conclusively  identir'y 
"bad"  data  since  it  is  extremely  difficult  te  distinguish  leasurenent  error 
f roi.i  grid  warpage).  -  gbgective  1. 


Inert  1  fled  "nan"  data  points  should  be  removeu  piendi  ng  further  aata 
g.itueriog  jod  i  nvout  .ga  ti  on  in  the  im.-ieiiate  vicinity  of  tne  points  in 
i  IPS',  ion.  -  Onjective  1. 


3.  "Bad"  aat 

a  points 

.lay  uopear  as  a  result  of  poor  initial  selection  of 

:  ‘i 

c. 0 13 1 a  i-  f  0 n  1 

pedu.nce. 

-  Objective  0. 

4.  A..'\:\P 

u'.  i  y  ■■ic 

-lA.*  1  s 

throughout  a  la.-je  range  of 

/a  1 

.■a>  “or  index  ' 

.,1^  'uq-rj.. 

t  ;.,)n  ind  vertical 

laps..',  bijt  these  progncjm  constants 

1 1 J 

.)e  '.ne  SdniO 

;n  the  AN 

/  A  \  \  .  si  i .  [  n  c  V  n  s 

stent  vertical  lapse  factor  can 

oro't.cu;  p.ositional  errors  greater  titan  51)0  ‘'eet.  Inccns  istiint  incex  of  refraction 
values  can  contribute  increnental  positional  errors  u^i  to  100  feet  jier 
staLior.  -  Objective  3. 

5.  The  oest  positional  acc  .racy  is  acnieved  when  tne  ei 
.10  A';/APN-101  matenes  the  modtl  used  in  XA-IP.  An  inconsist 
spiioroiii  mo'iels  for  tlie  input  oata  and  WAIIP  will  not  leasura 
j'os  I -Vial  uccijr  ]Cy.  -  Objective  3. 

3ri0jsm,_;  jn  j,.:,;r,jper  j^rth  nodel  'in  the  AN.'AilN-lOl  j 
:.jnt''ib,..i.i.:  uu  ti  104  fee.:  of  acditiondl  cos’tioncl  errjr.  1 
selectee  wnen  l.Cs/C  ^nouln  oe  use.i).  -  Oo^ective  3. 

.  The-  uAif’  i'’p..L  na  ,a  anc  Loren  ch.-.i'i  'iata  shoulci  be  relative  id  the  sa.ie 
system  Counli  uj-os.  -  Abjective-  3. 

0.  A  1)0  X  ,  n.ei:icil  .lie  prime  -lovee.^jj  area  with  tiie  urifor;,  data  point 
i-ins  ty  of  1  pr.'i'i:  i.;  •  n  na. it'cul  ..mle  cell  is  nighly  roconnienJoii.  ileducmg  I'lO 
<!  :'  1  .;e  .  y  ,  ;  lb  X  .  na,.;,i.:.il  ni  1  es  produced  a  co.iposite  positional  er-'or  of  7,'g 
rou.  ,  ./oile  ,1a  1  lit  ; .  ;i .  ■  ,  a  K  '  nautical  .niie  density  produced  only  a  114  foot 
i  error.  -  3;i^ect,ive  4. 

g.  .'.'ARi’  is  cnrreiitly  limited  to  400  input  data  points.  "here  is  no  physical 
I  .i;it,.it;on  ti  tne  :  n  nf  tne  f'CA. 
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-  so^s  .H  \var,)e.je  coei  r  ic  -or  sysce:n 

^)‘et,  3500-7500  feet,  /500-i:l501  feet,  and  above  12500 
ins : ,;r;ul  accuracy  will  Oe  aciiieved  with  4  independent 
i'lcted  v/ithin  tile  uppropi'i  ute  ..Ititude  ranges.  AN/AilN-lOl 
:  s  ..i.'.y  ie  satif’.ed  w’ t.n  r.upl  icate  or  i-edundant  .lata, 

’....ran  riav g.j i  ona  1  accuf'acgv  will  result.  'ost  results  indicate 
ei.jl  error  of  11)/  feet  throughu.it  i.he  OCA  from  0-15000  feet  wiien 
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using  4  indei’ondent  altitude  sets  of  warpage  coefficients.  When  only  one  set  of 
coefficients  was  used  (i.e  cither  the  IK,  5K,  lOK,  or  15K  data)  for  all  system 
altitudes,  the  resulting  average  error  ranged  froii  127  to  158  feet.  -  Objective  •>. 

11.  In  the  LORAN/LORAN  AN/ARN-101  navigational  node  the  system  may  experience 
instanteous  i)Osition  jumps  of  greater  than  1000  feet  when  transitioning  the  t'CA 
boundary.  In  higher  order  navigation  modes  Kalman  filtering  of  INS/LORAh  positijn 
will  prohibit  the  instanteous  jumps.  -  Objective  7. 

12.  When  a  saltwater  propagation  model  was  used  (i.e.  No  PROP  K  in  the 
AN/ARN-101),  the  average  positional  error  throughout  the  PCA  was  1261  feet.  The 
"Millington  method  of  deternining  a  single  constant  i-ipedance  value  reduced  r.lie 
overall  error  to  833  feet.  The  ideal  constant  value  (Rest  Ave)  further  reducmf 
tiio  error  to  523  feet.  When  utilising  the  full  warpage  polynomial  model  the 
overall  positional  error  vvas  117  feet.  -  Objective  3. 

13.  Software  prograni  modi f ications  were  i ncori)ordted  throughout  the  period  of 
this  test  activity.  The  SOW  (Attachaient  R)  lists  the  detailed  tasks  that  were 
accomplished  and  a  brief  highlight  of  the  major  jjrogram  changes  includes  the 

foil owi ng ; 

*  Deuailed  program  comments 

*  Improved  software  coding 

*  Extensiv(?  program  diagnostics 

*  Impro/ed  prediction  error  summaries 

*  New  plotting  and  data  smoothing  techniques 

*  Comi'rehensive  user  and  analyst  manuals 
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niiriiii,)  prep.ir.it,  ion  of  Lois  test  re))ort  an  improved  version  of  tlie  Warpaje 
fiic'f  r  ic  lent  Generation  Provjrain  (WARP)  was  nearing  completion.  As  a  result,  many 
of  the  (irovjraiii  recommerKlat  ions  that  would  have  been  described  in  this  section  will 
(le  i ncorjiorated  in  the  revised  program  and,  tiierefore,  will  not  be  discussed  here. 

duriiKj  tins  test,  effort  several  areas  iwre  note!  for  improving  ttie  AN/ARW-ll)l 
softwire.  It  is  reciKiimetuled  tiiat  future  efforts  lie  made  to  streamline  the  data 
rs.'gu  i  I'ei-.eiUs  of  Die  1.0RAN  ini  tial  i  j:at  ion  structure  within  the  AN/ARN-lfll.  The 
interaction  required  between  the  WARP  user  and  the  AN'/ARN-Ull  operator  is 
relatively  tiigh.  Unless  f.irtlier  vvorK  simplifies  data  inputs  to  the  AN/ARN-lOl, 
the  tr'.insfer  of  WARP  data  will  remain  subject  to  operator  error. 

It  IS  recotiiMiinded  lhat  the  'easibility  of  eliminating  the  i  nteriiedi  ate  step  of 
using  impedances  to  detenaine  Loran  warpage  be  i  nvest,  i  gated.  Preliminary  tests 
itiditMlr  (  ha(.  W'ARP  can  !ie  sii.iply  ,iiOi.lified  to  moilel  warjiage  in  time  units  instead 
.'f  impedance  with  no  loss  of  accuracy.  As  aiany  as  300-500  words  of  AN/ARN-101 
core  ..usaory  nay  be  saviivl  by  eliminating  the  impedaiks’  loot-up  tables.  A  small 
savings  in  AN/AR.N-ldl  progra-i  Oveculion  could  also  be  expected. 

■Inch,  analysis  roiiiams  to  be  tlone  pending  receipt  of  additional  operational 
lilt  a.  It.  is  iirobable  that  t  iie  input  data  density  can  be  expressed  as  a 
i,r  hemal  ical  limit  and  not  as  a  geographical  one.  ’'Ins  can  signif  icantly  reduce 
:  ir  data  collection  regut remeni-  It  is  also  probable  that  the  data  collection 
requ :  roiiiea;  s  as  a  funcl  i>)n  of  a.titude  can  tie  lore  clearly  delineated  if  the 
eflis.'ts  of  ‘errain  are  bid.ter  understood,  lailure  to  extend  the  analysis  will 
force  vno'y  Cvin.serva  1 1  vo  and  i'erhups  redund.int  data  collection. 

I'  IS  recoiiiniended  that  the  ilility  of  WARP  to  inlroduct'  "real  world"  warpage 
data  la  1  nr  ;  orr.n  grid  predict!  in  process  be  i  nvest  i  ga  ted.  At  present,  !.he 
ptMcrs  .  a;  collecting  "real  wiir  d"  l.oran  grid  data  and  subsequent  operation  of 
W\'tP  1..  I'l  iiKlepi'iident  eft’or't  I'rom  Loran  grid  prediction.  This  testing  has 
viol  0  r; :  that  WARP  lias  esser.t  ally  no  [irediction  capability  even  tliough  the 

uarpage  process  ’s  "s.nooln"  by  illing  in  nei  glnor ,  nj  data  for  void  area.  Rata 
o 't  fra,;  "tVieiutly"  territory  could  'oe  used  to  initiate  the  prediction 

,r\),  in!  1  "lint  r  iividly"  tern  ory  and  alsii  to  aiii  ;n  verifying  ttie  prediction 
ri",,,':  t  •  , ,,  -..•roisgly  ‘'nit  t  lat  an  eventual  ..larrying  of  the  WARP  program  and 

'  Mc  ic;,;  I  rogrj  1  is  'e.ilica  'lt’.  A  treaendous  savings  in  current 

coiip.'i'f  opera;  MO  tor  t  tu-  prOvI  ct  ion  pregt'o.'i  could  tie  acliieved.  In  addition,  tlie 
.ui’CiI:  ;r-iif  data  i.i  ,age.i;-”it.  sy.tem  would  In'  enfianced. 

Alt  oougl;  t  tie  current  Vto-sioi  of  WA'v'’  nas  been  tailored  to  aid  AW/ARN-101  Loran 
d.o  .1  rev,.n  rement  ,  inaninuiiig  oierat.ional  support  is  required.  It  is  felt  tiiat 
t  ne  opn-a;  lonai  coniands  will  n.'od  technical  assistance  in  establishing  prime 
civi'rir,e  nrea-i,  ilet  eiia;  in  ng  m'qiired  al  I  i '  adr  levels  tor  data  collection, 
rsi  I  'l'  IS)  tile  '.eventy  of  expected  loran  grid  waiqiage,  aid  cOfiputing  interaied i ate 
ronu  .ri'  impedence  v-iKu's  (tn  lieu  of  full  polyro.iii  al ) .  It  is  tiierefore 
imprra’  ve  Miat  a  mil  ra1  offic..'  be  est.ihl  i  shed  as  i  lie  focal  point  for  engineering 
'...ppii’d  !o  ,iid  t  le  liui'T'  ill  Pi.ita  g.it  tu’ri  ng ,  d.ata  r'Ovluclion,  and  analysis  of  WARP 
Oil!  p. Is.  il'ily  :>y  .'rovidi  i;  Slic  i  .i  grid  warpage  .lanager,  can  we  expect  to  actiieve 
'he  .pecitied  A.A.'Au.R-i  )1  pi'i'for .lanci'  accuracy  tor  navigation  and  woa)Hin  delivery. 
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Abstract 


Differences  between  LORAN  time  differences  (TDs)  measured  with  a 
receiver  and  TDs  calculated  using  secondary  phase  corrections  from  a 
homogeneous  LORAN  conductivity  model,  can  be  interpreted  as  i rregu- 
larities  in  the  actual  hyperbolic  lines  of  position  (LOP).  This 
condition  is  referred  to  as  warpage.  LORAN  warpage  occurs  for  two 
reasons.  First,  warpage  occurs  whenever  the  propagation  medium  is 
not  uniform.  Second,  warpage  occurs  when  LORAN  signals  travel  over 
two  or  more  different  propagation  mediums  i.e.  from  land  to  sea  and 
back  to  land.  The  presence  of  Loran  warpage  affects  the  normal 
coordinate  conversion  relationship  between  the  regular  LORAN 
hyperbolic  grid  and  the  corresponding  geodetic  (latitude,  longir  •;) 
grid  commonly  used  for  navigation. 

The  effect  of  Loran  warpage  is  corrected  by  the  AN/A  '1-101  L.RAN 
system  in  a  two  part  algorithm.  The  first  part  operates  off-line  and 
uses  a  paired  data  base  of  measured  LORAN  and  geodetic  coordinates  to 
generate  a  set  of  15  coefficients  for  each  LORAN  station.  These 
coefficients  are  used  in  the  second  part  of  the  algorithm  ear^ 
coordinate  conversion  cycle.  The  warpage  is  corrected  during  '■oordiratt. 
conversion  by  using  the  coefficients  tr  estimate  an  effective  wave 
impedance  for  the  secondary  phase  correction  of  each  LORAN  time  of 
arrival.  This  paper  describes  the  considerations  and  methodology 
needed  to  obtain  a  useable  data  base  that  meets  the  requirements  for 
coordinate  conversion  accuracy.  The  recent  calibration  of  the 
Southeast  U.S.  LORAN-C  chain  at  Eglin  APB,  Florida  is  used  as  an  example 
of  the  procedure. 
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Introduction 


For  many  app 1 i cat  ions ,  the  usefulness  of  Loran  depends  on  conversion 
from  hyperbolic  navigation  to  geodetic  navigation.  The  process  of  coordinate 
converting  Loran  time  differences  to  latitude  longitude  can  introduce  a 
significant  amount  of  positional  error.  This  error  is  caused  by  the  simple 
fact  that  the  Loran  grid  is  not  directly  related  to  the  geodetic  grid. 
Non-cancelling  signal  propagation  errors  cause  a  shift  or  bend  in  the 
smooth  regular  hyperbolic  lines  of  position.  This  is  called  Loran  warpage 
and  is  a  function  of  the  propagation  path  or  medium  the  Loran  signal  crosses. 

Loran  warpage  can  be  corrected  during  coordinate  conversion  by  using 
an  appropriate  model  for  the  secondary  phase  correction  term.  Most  coordinate 
conversion  algorithms  assume  a  regular  homogeneous  Loran  grid  and  model 
a  single  type  of  propagation  medium  for  the  entire  coverage  area.  This 
assumption  is  not  valid  for  accurate  navigation  over  land  where  more  than 
one  type  of  propagation  medium  is  involved  with  a  Loran  signal  path.  For 
very  accurate  coordinate  conversion  over  land  areas,  the  influence  of  each 
propagation  medium  must  be  considered.  In  areas  of  severe  warpage,  such 
as  mountanous,  this  requires  a  detailed  piecewise  addition  of  all  contributing 
factors  to  arrive  at  an  effective  delay  or  impedance  value.  This  process 
must  be  repeated  for  all  three  Loran  stations  (M,  A,  B)  in  the  triad. 

For  regions  wriere  less  coordinate  conversion  accuracy  is  required,  average 
impedance  values  can  be  selected  for  each  Loran  station.  These  values  are 
selected  based  on  the  average  propagation  medium  for  that  signal  path. 

Tne  U.  S.  Air  Force  has  demonstrated  very  accurate  coordinate  conversion 
over  areas  of  severe  Loran  warpage.  This  procedure  uses  a  two  part  algorithm 
for  determining  the  secondary  phase  correction*.  This  algorithm  is 
implemented  in  the  AN/ARN-iOl  Digital  Avionics  System.  The  first  part  of 
the  algorithm  uses  a  paired  data  base  of  measured  Loran  time  differences 
and  geodetic  latitude  longitude  coordinates  to  generate  a  set  of  coef r i c i en ts . 
The  coefficients  are  then  used  in  the  second  part  of  the  algorithm  each 
coordinate  conversion  cycle  to  compute  the  secondary  phase  correction  term. 
Thus,  the  accuracy  of  the  final  coordinate  conversion  depends  on  the  selection 
and  correspondi ng  accuracy  of  the  measured  data  base.  This  paper  presents 
tne  gi;nerjl  procedure  for  obtaining  an  accurate  data  base  for  the  AN/ARN-101 
coordinatt'  conversion  algorithm. 

General  Data  Requirements 

The  LORAN  coordinate  conversion  data  base  consists  of  LORAN  time 
differences  (TDs)  and  corresponding  geodetic  positions  (Lati tude/Longi lude) . 
Measured  (observed)  data  is  desired  but  interpolated,  derived,  or  calculated 
data  can  be  used  if  it  meets  the  accuracy  requirements  of  the  system. 

Tie  warpage  coefficient  generation  program  requires  that  five  geographic 
areas  be  defined  as  shown  in  Figure  1.  The  prime  coverage  area  is  centered 
over  the  desired  operational  region.  The  size  of  the  prime  area  is  variable 
but  it  Is  normally  limited  to  IOC  nautical  miles  by  100  nautical  miles  or 
smaller.  Within  the  central  prime  area,  the  warpage  coefficient  generation 
program  uses  the  LORAN  data  base  in  a  regression  model  to  calculate 
coefficients.  These  coefficients  are  then  used  in  an  interpolation  equation 
which  calculates  the  effective  wave  impedances  as  a  function  of  position 
in  the  prime  area. 
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I  Arei  5  4 

I  Northwest  Northeast 


Figure  1.  AN/ARN-101  Warpage  Correction  Areas 


The  un i f orm  distribution  of  data  points  witfin  the  princ  area  is 
very  important.  Ideally,  when  the  prime  area  is  divided  into  square 
cells  which  are  5  nautical  miles  by  5  nautical  miles  in  size,  each 
cell  should  contain  at  least  one  data  point.  This  distribution  insures 
that  the  prime  area  will  be  accurately  modeled.  This  is  a  minimum 
density  requirement.  More  data  points  per  cell  in  the  prime  area  will 
generally  result  in  a  more  accurate  mode  .  Locations  outside  the 
designated  prime  area  are  used  to  calcul  ite  an  average  wave  impedance 
for  the  four  areas  outside  the  prime  area.  Usually  data  points  at  a 
distance  of  10  to  20  nautical  miles  from  the  edges  of  the  prime  area 
are  sufficient  for  a  good  effective  impedance  model  for  the  outlying 
areas.  More  data  points  may  be  required  if  severe  warpage  is  encountered. 

The  AN/ARN-101  LORAN  warpage  model  was  designed  for  aircraft  and 
requires  LORAN  warpage  coefficients  at  four  altitude  levels  in  the  prime 
area.  Each  altitude  level  must  have  a  separate  data  base  according  to 
the  fol lowing  table: 


Leve  1 

Alt! tude  Range  for 
Data  Point 

Ground 

0  to  2500  feet 

5000  feet 

2500  to  7500  Feet 

10000  feet 

7500  to  12500  feet 

15000  feet 

12500  to  17500  feet 

Data  points  for  geodetic  locations  at  other  than  ground  level,  are 
not  required  to  be  the  same  geodetic  locations  as  the  ground  level 
point.  However,  the  requirement  for  a  m-nimum  of  one  data  point  per 
cell  does  apply.  For  other  apo  1  i  cat  ■ '^'s  such  as  vehicle  monitoring, 
only  ground  level  data  is  recuired. 
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Data  Base  Requirements 


This  section  presents  specific  requirements  to  develop  a  LORAN 
warpage  data  base.  The  Southeast  U.S.  LORAN-C  chain  at  Eglin  AFB 
Florida  is  used  as  an  example,  and  the  following  procedures  will  be 
covered . 


•  Selection  o^  the  Prime  Area 

•  Selection  of  Secondary  Stations 

•  Estimate  of  LORAN  Warpage 

•  Data  Collection  Requirements 

•  Generation  of  Warpage  Coefficients 

Selection  of  tne  Prime  Area 


The  Eglin  AFB  prime  coverage  area  used  with  the  old  East  Coast 
LORAJ-l-C  chain  was  a  42  by  72  nautical  mile  rectangle.  The  size  of 
tnis  prime  area  requires  a  minimum  of  126  data  points  evenly  distributed 
over  the  entire  area.  The  previous  data  base  provided  126  data  points, 
but  they  were  not  evenly  distr  buted.  These  data  points  were  concen¬ 
trated  over  the  west  and  east  :and  ranges  of  Eglin  AFB  with  a  scattering 
of  points  in  the  other  areas.  This  uneven  distribution  of  data  points 
caused  problems  with  coordinate  conversion  at  some  edges  of  the  prime 
area.  From  this,  it  Is  concluded  that  it  is  better  to  reduce  the  prime 
area  to  a  manageable  size  and  collect  data  points  that  are  evenly  dis¬ 
tributed.  Figure  2  is  a  map  of  the  Eglin  AFB  land  ranges  showing 
approximately  a  35  oy  52  nautical  mile  rectangle  located  by  the  North¬ 
east  and  Southwest  corners  shown  below. 

N£  30°  50.00'N 
85°  55.00'W 

SW  30°  15.00'N 
86°  55.00'W 

This  represents  a  reasonable  sized  prime  area  for  the  Southeast 
U.S.  LCRAN-C  chain  at  Egnp  AFB,  Florida.  The  map  is  delineated  by  a 
g r  i  0  every  5  minutes  of  angle.  Thus,  if  each  square  of  grid  contains 
at  leasr  I  calibration  point  the  minimum  distribution  requirements 
have  been  satisfied. 

Selection  of  Seconaary  Stations 

The  Southeast  U.S.  LORAN-C  chain  offers  different  combinations  of 
stave  stations  which  can  provide  coverage  for  Eglin  AFB.  The  complete 
data  for  this  chain  is  contained  in  Table  1. 
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TABLE  I 

SOUTHEAST  U.S.  LORAN-C  CHAIN  -  RATE  7980(SL2) 


Station 

Latitude  & 
Longitude 

Station 

Function 

Coding 

Delay  & 

Baseline 

Length 

Radiated 

Peak 

Power 

Malone, 

Florida 

30-59-38. 74N 
8S-10-09.3(M 

Master 

-  -  - 

1.0  MU 

Grangevllle, 

Louisiana 

30-43-33. 02N 
90-49-43. eOU 

W 

Secondary 

ll.OOOiiS 
1809. 54us 

1.0 

Raymondville, 

Texas 

26-31-55. Old 
95-50-00. 09W 

X 

Secondary 

23,000us 

,4443.38iiS 

'400  kU 

Jupiter, 

Florida 

27-01- 58. 4«1 
80-06-53. 52W 

Y 

Secondary 

43,000us 

2201 .88tis 

300  kU 

Carol Ina  Beach, 

R.  Carolina 

34-03-46. 04N 
77-54-46. 76N 

_ 

Z 

Secondary 

59.000us 

2542.74iiS 

700  kU 

Computer  generated  mpas  showing  the  Geometric  Dll  lotion  of  Precision 
(GOOP)  for  all  master-slave  combinations  of  the  Southeast  LORAN-C  chain 
were  computed  to  help  select  the  best  secondary  stations.  Analysis  of 
this  data  indicates  that  the  Raymondville  slave  (X)  generally  provides 
a  lower  GOOP  over  Eglin  AFB.  But  its  baseline  goes  directly  through 
the  selected  prime  area.  This  will  result  in  rapid  changes  In  GOOP  over 
the  Northeast  corner  of  the  Eglin  prime  area.  For  this  reason  the 
Ma lone-Grangev i 1 le-Jupi ter  triad  (M-W-Y)  is  selected  to  provide  the  best 
overall  coverage  for  Eglin  AFB.  Figure  3  shows  the  location  of  Eglin 
ArS  in  relationship  to  the  Southeast  U.S.  LORAN-C  chain  (7980)  and  the 
ole  East  Coast  LORAN-C  chain  (9930). 

Est  mate  of  LORAN  Warpage 

Rough  estimates  of  the  effective  impedance  for  each  station  are 
needed  to  determine  if  warpage  in  the  prime  area  is  mild  or  severe. 

These  estimates  are  also  initijlly  used  in  the  AN/ARN-101  for  all  areas 
and  at  all  altitudes  until  a  sufficient  data  base  is  available  and 
warpage  coefficients  are  computed. 
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These  initial  estimates  of  effective  impedance  can  be  made  using  the 
following  procedures; 


(1)  Using  a  map  that  covers  the  prime  area  and  the  three 

transmitters  (Master,  Slave  A,  Slave  B)  ,  select  a  point 
approximately  in  the  center  of  the  prime  area. 


(2)  Draw  the  line  from  the  selected  point  to  the  Master 
station  and  estimate  the  fraction  (x  )  of  this  path 
that  is  seawater,  the  fraction  (x  ^  land^ 

this  path  that  is  nornal  land,  anS^ERI  friction 
(^rough/dry  land)  of  this  path  that  is  rough/dry  land. 
Note  that  +  x^ormal  land  "*■  Xpough/dry  land  1- 

The  estimated  impedance  (i(^estimate)  is  then  calculated  by 
eva 1 uat i ng  ; 


(3) 


Repeat  step 


(0.001055)  (x  )  + 
sea 

(0.03)  (x  ,  ,  ^) 

normal  land 

(.05)  ^^rough/dry  land 
2  for  Slave  A. 


) 


(k)  Repeat  step  2  for  Slave  B. 


The  three  estimates  can  be  used  in  the  AN/ARN-101.  This  is  accomplished 
for  the  prime  area  by  setting  the  constant  term  a  to  the  estimate 
of  eftectivo  impedance  and  setting  the  other  1^  warpage  coefficients  for 
each  transmitter  to  zero  (a^  through  .  This  must  be  done  at  each 
of  the  four  altitude  levels.  This  procedure  provided  the  following 
estimates  of  effective  impedance; 


Mas  te  I'  -  Malone,  Florida 
lOOe  normal  land 
2l_  =  (0.03  *  1)  =  0.03 

'"m 

Secondary  W  -  Orangeville,  Louisiana 
■}l  sea 

95^  normal  land 

=  (0.001055  •  .05)  +  (0.03  -95) 

w 

=  0.0285 


Secondary  Y  -  Jupiter,  Florida 
55V  sea 

m5V  normal  land 

«  (0.001055  *  .55)  +  (0.03  •'^S) 

y 

=  0.0\k^] 
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The  data  in  Tab'e  II  was  tai-.ulatoci  ti-  ‘-sti'a'e  the  expected  positional 
accuracy  of  the  AN/ARN-10’  .('O'd  1  n.ue  conv<  Me’'  ‘^actor^’'',  the  oath  of 
the  t  i  fT)e  difference  line  of  position  ^  cOP  ^  i  re-,’ 1  ,in )  ,  and  relative  ranee  and 
bearinq  to  each  station.  This  data  shows  rxco'lent  ueometric  coveraqe 
for  the  Egtin  AFS  area.  The  i ns tant aneous  (plane'"')  path  of  a  "^D  LOP  was 
calculated  from  the  azimuths  to  the  .DRAN  stations  by  the  following 
eauat I  on : 


LOPs  = 

Tan"' 

cos  !//  c 

“ 

- 

s  ;  n  'p 

whe  re  ; 

Bear i nq 

to 

•'iaster  t  fat  ion 

II 

Bearinq 

secondary  s  t a  t i on 

Figure  3  shows  the  TD  Lines  of  Position  from  Table  II  overlaid  on  a 
map  of  Eglin  AFB.  Based  on  these  ca'culatlons  it  is  concluded  that  the 
LORAN  coverage  from  the  Southeast  U.S.  L0R>'^-C  chain  will  be  uniform  with 
TD  crossing  angles  of  approxi mate  1 y  72  to  77  degrees. 

Data  Collection  Require me n t s 

The  priTie  area  's  usually  sufficiently  large  that  the  data  base  wi!' 
be  collected  over  a  long  period  o'  time.  Th  measurements  taken  fo’’  the  data 
base  must  be  compensated  ‘^or  weathe’’,  diurral,  and  seasonal  effects  in  order 
to  be  consistent  and  accurate.  Ground  monitors  are  needed  to  accumulate 
a  history  of  TD  data  and  establish  standard  IDs  associated  with  the  monitor. 
TDs  collected  in  the  field  are  then  corrected  hy  an  amount  determined  from 
the  following  equation; 


*Q  factor  is  the  absolute  value  of  the  determinate  of  the  coordinate  conversion 
gradient  matrix.  Absoivit”  va'ues  qr^at-”-  n.!  are  required  for  conversion 
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TABLE  II 


TD  =  TD  +  TD.  -  TD  , 
meas  in^t  std 


where:  TD  =  Raw  measured  TD  from  the  field 
meas 

TD.  =  instantaneous  TO  of  the  Monitor 
I  nst 

TD  ^  ,  =  Standard  TD  of  the  Monitor, 

s  td 


The  data  base  may  be  collected  by  ground  mobile  receivers  at  knov/n 
benchmaks,  geodetic  survey  point:.,  or  by  an  aircraft  equipped  with  LORAN 
and  other  equipment  for  independent  ground  position  location.  When  aircraft 
are  used  to  collect  the  data  bast;,  the  instantaneous  aircraft  locations 
should  be  accurate  to  150  feet  (iCr)  or  better.  Measurements  with  position 
uncertainties  of  up  to  300  feet  ( Itr)  may  be  used,  but  only  if  more  accurate 
measurements  cannot  be  obtained.  When  positional  measurements  cannot  be 
made  to  accuracies  better  than  300  feet  ( Icr)  ,  then  the  data  point  should 
not  be  used.  If  this  condition  is  prevalent  over  the  prime  area,  then  the 
prime  area  cannot  be  adequately  modeled  and  “average"  values  of  effective 
impedance  will  usually  provide  the  same  statistical  accuracy. 

Generation  of  Warpage  Coefficients 

Collection  of  the  data  base  is  only  the  first  step  in  the  procedure 
to  correct  for  the  effects  of  LORAN  warpage.  The  data  base  must  be 
processed  by  the  warpage  coefficient  generation  program  to  caluclate  a 
set  of  fifteen  coefficients  for  each  LORAN  secondary  station.  This  program 
requires  the  following  data  as  an  input; 

•The  Northeast  and  Southwest  corners  of  all  5  correction  areas 

•LORAN  station  locations 

•Slave  emission  delays 

•Spheriod  Model 

•Transmission  ranges 

•Program  constants 

-  Effective  impedance  of  master  station  (estimated  average) 

-  Atmospheric  vertical  lapse  factor  (0.85) 

-  Atmospheric  index  of  refraction  (1.000338) 

•Data  base  for  the  area  of  interest. 

Currently,  the  warpage  coefficient  generation  program  is  in  IBM  card 
format  and  the  above  data  is  entered  on  IBM  punch  cards.  The  format  for 
this  data  is  specified  in  the  program  users  manual*. 
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The  input  data  base  is  used  by  the  warpage  coefficient  generation 
program  to  model  the  warpage  and  estimate  the  effective  impedance  for 
coordinate  conversion  of  each  data  point.  Once  the  effective  impedances 
are  determined  to  the  accuracy  of  the  data  base,  warpage  coefficients  are 
generated.  The  current  version  of  the  program  uses  only  one  coefficient 
to  mode)  the  master  station.  The  estimated  value  of  master  impedance  used 
for  input  is  also  the  coefficient  for  the  master  station.  For  each  secondary 
station,  one  of  the  15  coefficients  is  the  average  effective  impedance 
for  that  station  while  the  remaining  l4  coefficients  are  constants  for  a 
4th  order  polynomial.  This  polynomial  is  used  to  adjust  the  average 
effective  impedance  as  a  function  of  pos i t  on  within  the  prime  area.  With 
this  method,  the  coefficients  can  effectively  model  the  warpage  and  provide 
the  wide  range  of  impedances  needed  for  correction.  This  process  is 
implemented  in  the  second  part  of  the  AN/ARN-lOl  coordinate  conversion 
al gor i thm. 

Summa  ry 

For  many  applications  LORAN  must  be  coordinate  converted  to  latitude 
and  longitude.  Changes  in  the  propagation  medium  over  land  signal  paths 
warps  the  smooth  hyperbolic  lines  o^  position  and  at^fects  the  accuracy 
of  LORAN  coordinate  conversion.  The  effects  of  warpage  can  be  correced 
by  proper  selection  of  the  secondary  phase  correction  term  but  the 
selection  is  complicated  by  the  severity  of  the  warpage.  Mild  v;arpage 
can  be  approximated  by  a  regular  smooth  LORAN  grid  ana  a  single  effective 
impedance  correction  model.  Correction  for  severe  warpage  depends  on 
location  and  requires  an  individual  correction  for  each  station  and  each 
location.  Precision  coordinate  conversion  in  areas  of  mild  and  severe 
warpage  requires  an  accurate  data  base  of  known  geodetic  positions  ang 
measured  TDs.  This  data  base  is  then  used  to  model  the  LORAN  warpage.  !n 
the  AN/ARN-101  algorithm,  warpage  coefficients  are  developed  in  a  warpage 
coefficient  generation  program.  These  coefficients  transfer  the  waroage 
model  to  the  second  part  of  the  algorithm  which  selects  the  proper  secondary 
phase  correction  for  coordinate  conversion.  This  algorithm  has  demonstrated 
the  ability  to  accurately  mode!  both  mild  and  severe  warpage  as  well  as 
the  rapidly  changing  warpage  encountered  with  multiple  1 and-sea- 1  and 
i nterfaces . 
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A^'PENDIX  ;i 


STATF:.VEf,'T  OF  WORK 
IlPDATF  OF  DOCUnEi^TATION  OF 

•■rpa;'f.  coefficient  oener/jion  program 


1.^  nJTRnDjjrjin': 

Tho  thp  "arpaoe  Coefficient  Feneration  Program  (WCGP) 

tn  tr--’n':frrr  LO-A.i  nrid  calibration  data  (Loran  time  difference 
:■  at  kro'.-.n  oeccJetic  nrsitions)  into  ARN-101  warpage  c(  pffi- 

-^unction  is  necessarv  to  define  the  APiv-lOl  Loran  warpage 
' .el  tr  prcvii  e  nocessarv  navigation  ard  wcanon  delivery  accuracy. 

;..g  FCOP^ 


Tnis  effort  covers  tiie  update,  modification,  documentation,  and 
checkout  of  all  '■'rccrai'  functions  witiiin  the  U'arrage  Coefficient  Genera¬ 
tion  ^roorar^i.  All  non-renuired  program  logic  will  be  eliminated,  necessary 
Icoic  strca'.nl  ineu,  reouirod  capabilities  added,  nrogram  variables 
stan^ar.lized,  aoc  user  and  analyst  manuals  prepared  providing  complete 
and  uo-to-datc  dccur.;entation  of  tlie  program  capabilitv.  All  options 
retaiiied  witniii  the  updated  program  will  be  executed  by  means  of  a  sample 
nrcl.’l oil),  results  .tssessed,  and  if  required,  further  modifications  made 
tc  prrducc-  desired'  results. 

_0  dE‘,LhAL  k 

"he  a r Coefficient  Generation  t>roqrdri  vas  developed  by  Lear 
ieei'T,  Inc.  (i-Si;  as  a  support  program  for  the  AN/ARN'-lOl  Loran  function. 
It  i.as  provided  to  t.he  Armament  Division  (AD)  for  support  of  the  AN/APi;-101 
;r..';g/c.,:£  bT loj," E  test  procran..  Delivery  consisted  of  a  card  deck, 
so'joc^  listinn,  and  User's  ’-lanual  (LSI  YVlOOd).  The  program  resided  at 
until  dnverb'^n  l}7r.  Thn  orieiaal  LFI  oroqram  was  then  provided  to  HQ 
,d ,  i.anS'tcT  A"3,  ‘’A  where  it  was  made  operational.  ESD  Operating 
Locatinn-''F  at  "el  in  AFF,  Florica  3?oA2,  received  the  nrogram  in  1979  for 
)se  vdf;  t'’e  '''^ii  A'T-  calilTation  of  the  new  Southeast  USA  Loran  C  Chain. 

T^e  ,'ar'->->’e  Cc'^fficient  Generation  '’roqram  is  currently  operational, 
tut  its  opera  I' '■")  is  complex  and  reouires  evaluation  by  an  engineer  we1l- 
'.'ersed  ip  Loran.  The  current  imrlementetion  includes  subroutines  taken 
:1rrr;,lv  fnc’"  ptjirr  functional  oroqrams  which  results  in  unnecessary 
.jnlic,K,';r  i:f  mutines.  calculations,  and  extensive  arrays.  The  program 
's  :  u'.h  ur  -'r  ani.  takes  lorge''  to  execute  than  is  necessary  for  the 
‘■nyr.tiin  .UK.  eviluatiop  cf  v/d>  r  -.''ficients.  In  addition,  t'u* 
'n-nrai'  tnal/sis  ••  ruin'.'  for  su.tpssful  utilization  is  difficult  to  per- 
forn  do’-UMentation  available. 


I 


J 


T  'S:./Tcai;._icAL 

'  ,.^  c'  rtr.ictcr  .*1^11  nerfonn  tut-  specified  herein.  All  tasks 

s  kmi  ..I.  i  or''.;! t  ,0(1  luj  the  (  ericd  of  time  specified  in  the  schedule. 

Pi  ’c  i", '/ici  for  ossianed  tasks  shall  be  as  directed  by  the  government. 

'’.1  ‘vn'i/ize  arcj  siiinlify  input/output  data  reouirements  necessary  to 

i.t^-  the  prcnron’. 

V!.'  ’^rovid^  I  '^inition  of  all  input  parameters  with  format  specifications 

*:  ontT'v  si^nsitive  inout  parameters  and  provide  nominal  values. 

tc  ,'atd  point  input  onto  one  card  by  eliminating  the 
.  ■>;;  a>:..  roformatinq  the  remainino  data. 

.1.'.  I'.pfon  .'t  out-iut  to  duplicate  input  file  fonvidt. 

rovid'  ’('i’c  to  identify  and  eliminate  data  points  outside  the 

;  c  iC'  . 

•  •  1..  tra.  c' ."‘di nation  and  approval  of  the  government,  delete  all 

'  rrcy  jdi;  :  . ‘.,  ,1  it.'d  Icnic  and  streamline  required  logic  to  provide 

•^r.  ♦'d  Ic  .r  ;'  uo-  tc-date  programming  techniques.  All  code  must  be 

rill r  par  ANSI  STC  X  3.10  -  1978. 

.  ...i  ct1l  creation  process  to  accomodate  residual  regions 

it,  lir  i  ni'i:  -’  crver.Kio  area. 

.1...  ;  land. if  11  variable  names  and  reorganize  COMMON  to  achieve 

1  ;  y.  1  ir  routines  tc  provide  easier  flow  of  program  logic. 

K  .1..^  '  'i  -lion  of  program  variables,  variable  names,  units  of 
.r.  ,Tt.  a:  '  APoctfed  rungts. 

tifips  to  consistently  maintain  numerical  significance. 

. .  ’  ■  -  '  ■  '  'rpipctaLle  spheroid  models  and  allow  for  user 

'•  ■  ‘  ei.  neference  Specification  CPinol-004  of  Contract 

■  i  .V):  dcr.iputer  Proerar:  Development  Specification 

'  r  ^  '  ri,  '  :  i',pr  of  Pr-rlc  Dioital  Modular  Avionics  System. 

' :r  '  li  to  input  transmitter  elevation  and  use  it  in 


■  ;  'oritliin  to  coii'pute  all  areas  boundaries  given  only 

I  ri;  '  ,  '  I  ic  nrii’e  area . 
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■I.!.:'  i’adi^v  lU'-trix  conversion  ainorithr  to  convert  time  difference 
;7;i)  errors  to  units  cf  linear  measure. 

I'  rlo'  oi.t  liistoorar  niottinn  canability  for  selected  program 

'■  La . 

''.l.i?  Imnlenet::  additional  overlav  to  include  selectable  nlots  of: 

a.  Da*,-  'Hi'-ts  \'ithin  the  nrime  coverage  area  showing  grid 
■•■It'  C’ n  t'oundr’rics  and  cell  numbers. 

t'.  Co'i'-our  -'lots  of  calculated  and  estimated  impedances  for 
■:-':Cn  oaster/slavc  cor'bi nation . 

t.  Cor’lojr  nlot  of  imredance  error  measurements. 

■  .1.11  Add  to  tie  nroqraia  listina  comments  to  fully  describe  program 

ii  f:ic  and  cco.-^utaL n  ns  being  performed,  (nclude  description  of  program 
.otnii’v'’  assji  rt  ions  ,  and  mathematical  modelinn  techninues  employed). 

•  r.  ilc  'inc:  cv  oil  tior.  of  Task  4.1,  the  contractor  shall  conduct  a 
c  ictoit  of  all  ’^ro'rapi  controls,  options,  and  input/output  formats. 

■  L, . '  In  coorj  r, alie  n  with  and  with  the  approval  of  the  government, 

c  c, ''trace,  r  coall  prepare  saii:ola  problems  which  will  exercise  every 
■'r'onr.r  flomenl  'hcr.'bv  demonstrating,  to  the  satisfaction  of  the  govern- 
p'lt.  thu  into  ritv  i  '’  the  completed  program.  Test  data  will  be  supplied 
the  aj'v.  r': 'I '  "  re!,  later  than  30  days  prior  to  test  date. 

'■.1.2.'  roiicwir  ;  the  novernment ‘s  reviev;  of  the  results  produced  in  Task 
'.2.1,  the  coidr.iclor  shall  conduct  further  modification  work,  as  required, 
t  c  '•■'i'ct  a  !■'  r''i.a'nina  nroaram  deficiencies  as  identified  bv  the 
•’'•■vr'-n!  •-'P*-.  “(  ci'^ic.itions  may  take  the  form  of  additions  and/or  chanoes 

npoTp;"  lO'iir  ard  i I'put/cutput  formats. 

•  .  ■'■'0  cortractor  shall  conduct  a  documentation  effort  consistino  of 

r-'Plac'' "C'Pt  of  t 'e  nrpsent  L'arpane  Coefficient  Gercration  Program  documen- 
*iti'’r>  (L'^I  ‘lav  lO/C)  wit'""  the  followino: 

.  I  Jsm  's  daruicil  -  This  manual  shall: 

‘^La  .i-'  t-r-  rb-ipctive  of  the  computer  program,  its  potential 
!  f;l  ice  Lie  O'; .  an  J  nrovi'Je  concise  description  of  the  major  elements 
'  -'Pi*",  r  i|i  ;.i .  r  ri (i.e.,  overlays,  programs,  subroutines,  and 
■Id' ii  c  Li (  I'.s  } .  larr.itivp  discussions  and  illustrations  shall  be  used 

:  •,  T'[  t  '-jor  i  d  . 

'.L  osc"'  ■  .  '  '  :1  the  ■!  jitiLs  necessary  to  use  the  program  and 


f 


.  '^rocitv  ,-«,r,d  define  innut  narameters  reauired  for  each 

.  yt  ruticM  ant'  f  'C'ir’  location  and  format  on  the  input  data  cards. 

.  ’  Iliti'.trat-  data  utni-  >fiu'  -^or  all  major  options. 

.ale  the  output  options  available  and  define  the  output 

r  i.'it  tens  . 

define  output  items  givinq  dimensions  and  units  for  each 

/rt.  jpt'cn. 

.d  d.ir/i'  j-.ip-'les  of  eacii  cutout  option. 

niu'  Liot..'  and  describe  eacl',  output  plot  option. 

.'c  SLi.t  a!i  exantple  of  hov/  to  exercise  the  program  using  a  ! 

.  '.l.  '.i  desLtiie  tiic  prohlem  including  assumptions  and  constraints. 

( 

J.  !h,:lu.'e  a  tabular  listing  of  the  input  parameters  with  the 

'  rinjl  va’in.  -  ’.o  “0  used. 

.  '  i  1  lii'  T.  the  total  punched  data  card  setup. 

...  'i  ■)'.  St':  ".u:"..  •.!  -  This  manual  shall: 

‘  iectives  of  the  comouter  oroaram,  state  its  potential 
■  1  '  .  r'  -fnvide  .a  concise  description  of  the  major  elements 

■  ryl'.'St's  ‘'tanual . 

■  "  '(  tioii  v.'hicn  describes  the  application  considerations  : 

*  •  ■  '  ■  1  "  ,  •  '.■arnaqc  node! .  1 

'  t (  ’  >9  data  ccllection  renuirements . 

...  I  f  chp  ur.c jl ibratea  grid  operation. 

. :v  :  :  •  nionitorincj  procedures  to  be  used  to  determine 

■  1 1  :  .  ;  r  !  .  .  •  r;  i  . 

.  f  V  i  .  V  cior.  v.'i;ic;i  cescribes  the  mathematical  and  spheroid 
.  s.e  updaiec  computer  program. 

.  y  natneniatical  enuations  used  in  each  model. 

.  ..  .  ,  111  too  assumptions  used  in  eacli  model. 
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.  '.J  i'o'ii  '  all  svr.ibol?.  abbreviations  used  in  each  model  and 
•rovi  '  t' i  '.ritr  (>'  each  ccr:staTt  and  variable  used. 

'.L'.r.''  v''^'ordinate  s/'^tems  and  transformations  for 

’  I'l'fratiral  ,c' '  Is. 

•'  s'.-ctior  v;hicb  defines  the  basic  desion  ohilosoohy 
i'  ■  ;  '"^lomertation  scheme  aonlied  to  the  math  models. 

‘'roviJt  0  functional  block  diagram  indicating  the  major 
crTrrr’  ts  rf  ‘.'r  rroqram,  their  interrelationships  and  their  inter- 
•'■  ■-cs  .it:^  oth:  >'  coroononts .  Includes  overlays,  programs,  subroutines. 


‘  Prcvi  ;l  a  Ionic  dianram  that  depicts  the  structure  and  logic 

0-'  ^(.ic  cvurlav  routine  and  subroutine  contained  in  the  updated  computer 

!  1  i  i  ’  J'.i'  I  . 


•.  '.t'.'.i  Ircluue  a  alossary  of  symbols  and  terms  used  in  the  software 
cu'"‘  rtj  tion. 

‘•■l.t.s  ^rovi..t  )  section  v/hicb  corioletely  describes  the  procedures  for 
*  rc 1  i.r  cnc  i.ecu>  entino  the  prooram. 

“^rovi 'e  voerational  instructions  for  tl'.e  pronram.  Include 
I'O'-iit  ro-'i,i  1  irii  taticr.s ,  as  v.'ell  as  interpretations  of  all 

■'  s  s  ’ .  in  ' n.  lou  . 

l  os  r'le  Magnostics  and  explain  all  error  code  messages. 
Incii.  ■  listinn  p-r  comouter  program  in  source  lanouage. 
'■  i  1  0  all  source  card  formats  identi^vinn  their  puroose 

:  fu''C‘:’r:n. 


’.3,3  llcv'i'T  '_h'’  nevornmert's  review  cf  the  results  produced  in 
i  isV'  •  .;..l  a!',(l  -.'..r,  t^n  contractor  shall  conduct  further  modification 
'  ,  as  rcouir.  ',  tr  crrrect  erv  deficiencies  identified  bv  the  government. 

.  T'l:  ccntri'.ctnr  shaT’  install  and  exercise  the  computer  program  on  the 
i  C  tCv',  1  cor  puti  r  at  Lolin  APB,  FL.  This  task  will  be  completed  when  the 
c( hir'utc'r  ">ro, -rat  .  runninn  on  the  Foiin  conouter,  is  able  to  duplicate  the 
'  .tnut  results  '''tained  Uurine  t'c  I'ronram  checkout  phase.  Task  4.?,  above. 

1 '  u  ')  veriiLPt  arranue  for  tee  use  by  ttie  contractor  of  the  Eglin  AFB 
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